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INTRODUCTION TO COLLIDER PHYSICS 



MAXIM PERELSTEIN* 



Newman Laboratory of Elementary Particle Physics 
Cornell University, Ithaca, NY 14853, USA 
, E-mail: mp325@cornell.edu 

This is a set of four lectures presented at the Theoretical Advanced Study Institute (TASI-09) in June 2009. The 
' goal of the lectures is to introduce students to some of the basic ideas and tools required for theoretical analysis of 

collider data. Several examples of Standard Model processes at electron-positron and hadron colliders are considered 
to illustrate these ideas. In addition, a general strategy for formulating searches for physics beyond the Standard 
Model is outlined. The lectures conclude with a brief survey of recent, ongoing and future searches for the Higgs boson 
and supcrsymmetric particles. 
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qj ' 1. Introduction 

r- 1 | 

Our knowledge of the laws of physics in the sub-nuclear domain (at distance scales of about 10 -13 cm and 
smaller) is for the most part derived from analyzing the outcomes of high-energy collisions of elementary par- 
ticles. While the size and sophistication of each component of high-energy collision experiments have steadily 
grown, the basic experimental setup has remained unchanged since late 1960's. First, a particle accelerator 
uses a carefully designed combination of electric and magnetic fields to produce narrowly focused beams of 
energetic particles (typically electrons, protons, and their antiparticles). Then, two beams collide head-on, 
usually with equal and opposite momenta so that the center-of-mass frame of the colliding system coincides 
with the laboratory frame. a The region where collisions occur (the "interaction point" ) is surrounded by a 
set of particle detectors, which attempt to identify the particles coming out of the collision, and measure 
their energies and momenta. 

Since physics at subatomic distance scales is governed by laws of quantum mechanics, the outcome of 
each collision cannot, as a matter of principle, be known ahead of time; the best that any theory can do 
is to predict the probabilities of various possible outcomes. Modern collider experiments collect and analyze 
outcomes of huge number of collisions; the number of events with specified properties within the collected 
data set is proportional to the probability of such an event. Thus, the probability of a specific outcome 
of a collision (or a closely related quantity typically used in particle physics, the cross section) provide a 
natural bridge between theory and experiment. Narrowly speaking, one can say that a particle theorist's job 
is to infer laws of physics from experimentally measured cross sections. There is no known algorithm that 
can do this job. (This is lucky for us, since otherwise theorists could be replaced by computers!) There is, 
however, a well-developed formalism for predicting cross sections, given a Lagrangian quantum field theory. 
In practice, therefore, theoretical interpretation of collider data proceeds by picking a candidate theory, 
computing relevant cross sections within that theory, comparing with data, and (if the comparison does not 



a For about three decades before the advances in accelerator technology made it possible to steer two beams into a mid-air head- 
on collision, the principal technique was to accelerate a single beam and crash it into a stationary target. I will always consider 
a collider setup in these lectures, and all experimental data I will show arc from colliders; however, in terms of theoretical 
interpretation, the main subject of these lectures, the differences between the collider and fixed-target setups are rather trivial. 
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work) moving on to the next candidate theory. For the last 30 years, the default "leading candidate" theory 
has been the Standard Model (SM). As everyone knows, no statistically significant deviation from predictions 
of this theory has been observed so far, although at times mild inconsistencies with data have motivated 
theorists to try alternative candidates. It is, of course, equally well known to the TASI participants that 
strong theoretical reasons exist to expect that the SM hegemony will finally break down at the energy scales 
around a TeV, which will be explored experimentally for the first time by the Large Hadron Collider (LHC) 
in the next few years. Obtaining detailed quantitative predictions for the LHC experiments, both from the 
SM and from alternative candidate models, is a crucial task for theorists in the LHC era. The aim of these 
lectures is to introduce TASI students to some of the basic concepts and theoretical tools necessary to make 
such predictions. 



1.1. Definitions and Basics 

Consider a collision of two elementary particles, A and B, in the reference frame where the net momentum 
of the pair is zero. This frame is called the center- of-mass frame, or c.o.m. frame for short, since the center 
of mass of the system is at rest. In the case of e + e~ colliders, A and B are just the electron and the positron, 
and the c.o.m. frame coincides with the lab frame. For hadron colliders, A and B are partons (quarks or 
gluons), and the c.o.m. frame, which we will also call the "partem frame" in this case, is generally moving 
along the collision axis with respect to the lab frame. In either case, we will neglect the masses of A and B, 
since they are tiny compared to the energies we're interested in (of order 10 GeV and higher). By convention, 
we will choose the z axis to lie along the direction of the A momentum. The four-momenta of the colliding 
particles are 

p A = (E,0,0,+E), PB = (E,0,0,-E). (1) 

The total energy of the colliding system, the center- of-mass energy, is E cm = IE. We will also frequently use 
the Mandelstam variable s = (pa + Pb) 2 = E^ m . In the case of hadron colliders, the center-of-mass energy 
in a collision of two partons will be denoted by to distinguish it from the energy of the colliding hadron 
pair y/s. 

In particle colliders, collisions actually take place between beams containing large number of particles. 
If two beams collide head-on, the number of collisions leading to a final state with particular characteristics 
(type of particles, their momenta, etc.) should be proportional to the number of particles in each beam, Na 
and Nb, and inversely proportional to the beams' cross-sectional area A. The coefficient of proportionality 
is the scattering cross section for this particular final state: 

(Number of events) • A 
a= NaNb ' (2) 

If beams collide at a frequency / Hz, the rate R (number of events of a particular kind recorded per second) 
can be written as 

R = L ■ a , (3) 

where 

L = (4) 
A 

is the instantaneous luminosity. Simple as it looks, Eq. (3) is a fundamental cornerstone of collider physics, 
and it is worth examining it more closely. The rate R is measured directly by experimentalists. b The quantity 



b Actually, what is measured is R ■ £, where £ is the detector efficiency: the probability that an actual event with particular 
properties is identified as such by the detector. Efficiencies vary widely depending on the detector and the kind of process one 
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L (together with E cm ) contains all the information about the accelerator needed to analyze the experiment. 
Experimental collaborations carefully monitor and record L, as a function of time. The experimentally 
measured value of the cross section is inferred from Eq. (3). This value can then be compared with the 
theoretically expected cross section. 



Table 1. Recent and future energy-frontier particle colliders. (Parameters listed for the LHC and the ILC are design 
values.) 



Name 


Type 


(GeV) 


Lint (pb- 1 ) 


Years of 
operation 


Detectors 


Location 


LEP 


e+e~ 


91.2 (LEP-1) 


200 (LEP-1) 


1989-95 (LEP-1) 


ALEPH, OPAL, 


CERN 






130-209 (LEP-2) 


m 600 (LEP-2) 


1996-2000 (LEP-2) 


DELPHI, L3 




SLC 




91.2 


20 


1992-98 


SLD 


SLAC 


HERA 




320 


500 


1992-2007 


ZEUS, HI 


DESY 


Tevatron 


PP 


1800 (Run-I) 


160 (Run-I) 


1987-96 (Run-I) 


CDF, D0 


FNAL 






1960 (Run-II) 


6 K (Run-II, 06/09) 


2000-??? (Run-II) 






LHC 


PP 


14000 


10 K/yr ("low-L") 


2010? - 2013? 


ATLAS, CMS 


CERN 








100 K/yr ("high-L") 


2013?? - 2016??? 






ILC 




500-1000 


1 M??? 


??? 


??? 


??? 



Throughout the lectures, we will contrast theoretical predictions with data from recent and ongoing 
experiments at energy-frontier colliders. Table 1 shows the basic parameters of these colliders, along with 
the upcoming LHC and the proposed next-generation electron-positron collider, the International Linear 
Collider (ILC). C It is important to keep in mind that, for hadron colliders, the listed center-of-mass energy 
corresponds to the colliding (anti)protons. Since high-energy processes are initiated by partons, which only 
carry a fraction of the proton momentum, the energy scales that can be probed at a hadron collider arc 
substantially lower than this energy, typically by factors of 3 — 10 depending on the process. Electron-positron 
colliders, on the other hand, are able to explore many reactions at energy scales extending all the way to 
their nominal <Js. The luminosity values shown in the table are the integrated luminosities, L lnt = j Ldt over 
the lifetime of the experiment. The table also lists the detectors at each collider. Detector names coincide 
with the names of collaborations of physicists operating them, and are frequently used to refer to the data 
published by these collaborations. 

Computing and interpreting cross sections will be our main focus. It is clear from its definition, Eq. (2), 
that the cross section has cgs units of cm 2 . A unit typically used in experimental nuclear and particle 
physics is 1 barn = 10~ 24 cm 2 . In "theory units", c = h = 1, the natural unit for cross section is GeV -2 ; the 
conversion factor is 

1 bn = 2568 GeV -2 , 
1 GeV -2 = 3.894 • 1(T 4 bn. (5) 

To get a very rough estimate of cross sections expected in particle physics experiments, we can use dimensional 
analysis: away from thresholds and resonances, the only energy scale in a collision of two massless particles 



is considering. In addition, the measured rate typically includes events that do not actually have the requested properties, but 
are mis-identified due to detector imperfections. In these lectures, we will mostly not be concerned with such detector effects, 
except for an occasional brief comment. An interested reader is referred to Eva Halkidakis' lectures at this school. 
c For lack of time, I will not be able to discuss results from recent lower-energy, "luminosity-frontier" collider experiments, such 
as CLEO, BaBar, and Belle. 
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is E cm , and we should expect the (total) scattering cross section to behave roughly as 

1 



a 



E 2 



(6) 



A similar result (larger by n) is obtained by replacing the colliding particles with classical "billiard balls" of 
radius equal to their Compton wavelength A <~ 1/E, and taking their geometric cross section as an estimate. 
The geometric cross section also coincides with the upper bound on the total inelastic cross section (assuming 
s-wave scattering) from unitarity considerations. The cross sections for specific processes are typically lower, 
by an order of magnitude or more, than this bound: For example, the e + e~ — > Z cross section on resonance 
(•y/s = Mz) is about 40 nb, compared to <7 g0 om = v/M^ w 2500 nb. The decrease of cross sections with 
energy has an important implication for accelerator design: Colliders operating at higher center-of-mass 
energies must also have higher luminosity, adding to the technical challenges of expanding the high-energy 
frontier. This trend is clear in Table 1. 

The "master formula" for evaluating the cross section and kinematic distributions for a 2 — > N scattering 
process is 

da = Ys (6 ^ ak) • (2*) 4 * 4 (PA +Pb J>) • \M(p AlPB -+ fe}| 2 , (7) 

where M. is the invariant matrix element, a.k.a. scattering amplitude, and pi = (E^pJ are the 4-momenta 
of the final-state particles. Note that M contains all information specific for the process under consideration 
(such as coupling constant dependence, etc.), whereas all other ingredients are simply kinematic factors 
common for any 2 — > N process. While Eq. (7) is written in the center of mass frame of the colliding 
particles, it is in fact invariant under boosts parallel to the collision axis. This feature will be important 
when hadron collisions are considered. If the colliding beams are unpolarized, one needs to average the 
quantity \A4\ 2 over all possible initial-state polarizations. If the beams are polarized (this was the case at the 
SLC, and may be implemented at the ILC), an appropriately weighted average should be computed instead. 
In addition, if the final-state particles have spin, \M\ 2 should typically be summed over all possible spin 
states, since no collider detector is capable of detecting spins of individual particles. (Exception occurs when 
the final-state particles decay promptly, in which case the angular distribution of their decay products may 
carry information about their polarization state.) The appropriately averaged and/or summed Scattering 

2 

amplitude| 2 will be denoted by \M\ . 

The number of independent kinematic variables in a 2 — > N process is 3N — 4. In practice, the initial 
state is always symmetric under rotation around the collision axis, and no physical observable can depend 
on the overall azimuthal coordinate, leaving 37V — 5 physical variables. The simplest case, most commonly 
encountered in practice, is 2 — > 2 scattering. The only observable not constrained by energy and momentum 
conservation is the scattering angle 9, which by convention is defined as the angle between the 3-momenta 
of particles A and 1. The differential cross section is given by 



do { pyi tf y /s>m 1 +m 2 ; 



dcos9 l^o otherwise, 

where 



(8) 



| pi | _ 1 / (s-m\- m 2 ) 2 - 4m 2 m 2 2 ^ 

In the most common case of equal masses in the final state, mi = m 2 — m, this formula further simplifies 
(for y/s > 2m) to 



dcos9 32-7TS 



da 1 h-^Wt do) 
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Note that the square-root factor is simply the velocity of the final-state particles (in units of c) . The quantity 

\A4\ is often expressed in terms of the Mandelstam variables, Lorentz-invariant (scalar) bilinears of the 
4-momenta of incoming and outgoing particles. In the case of 2 — > 2 scattering, these are 



Pb) 



t = {pa- Pif 



P2) 



(11) 

The Mandelstam variables are 



They are not independent: it can be easily shown that s + 1 + u = m\ + 
related to the scattering angle: for example, if mi = m 2 = 0, we simply have 

t = -| (l-cosfl), u = -| (1 + cos(9) . (12) 

The main advantage of using Mandelstam variables comes in applications of crossing symmetry to relate 
processes such as, for example, electron-positron annihilation e + e~ — > 77 and Compton scattering e~7 — > 
e _ 7. They are also convenient for analyzing hadron collisions, being invariant under boosts connecting the 
parton and lab reference frames. 



2. Electron-Positron Collisions 

In this Lecture, we will study a few examples of reactions initiated by electron-positron collisions, and use 
them to illustrate some fundamental concepts and issues central to the field. Since colliding particles are 
elementary, e + e~ collisions are somewhat easier to analyze than collisions between hadrons, which will be 
discussed in Lecture 3. 



2.1. Muon Pair- Production 

We start with the process e + e~ — > At tree level, only two diagrams contribute, see Fig. 1, making 

it possibly the simplest 2 — > 2 reaction in the SM - the "hello world" example of collider physics. Most 
(probably all) TASI students would have calculated the cross section of this reaction in their Quantum Field 
Theory (QFT) classes, probably using four-component (Dirac) notation and trace technology to perform spin 
sums. An alternative is to use two-component (Weyl) fermions, and to evaluate the scattering amplitudes for 
particles in definite helicity eigenstates. This method provides more insight into the physics of the process, 
and becomes especially valuable when weak interactions are considered. Let us outline the calculation. 




Fig. 1. Leading-order (tree-level) Feynman diagrams contributing to the process e + e — > fJ. + fi 



2.1.1. Muon Pair-Production in QED 

To begin with, let us only consider the diagram with the virtual photon exchange, ignoring the Z. (This 
would be a good approximation at low energies, y/s <C Mz.) In two-component notation, electrons and 
positrons are described by a pair of two-component (Weyl) spinor fields, ej, and en. The subscript L/R 
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denotes the field's chirality, which determines its transformation properties under the Lorentz grou p. Th e 
QED Lagrangian in this notation is 



C 



ieU"Df;e L 



m e 



e R e L + e\e R 



Here = = <9 M — ieA^ are covariant derivatives, and 

a = (1,5), a = (1, -5) . 



(13) 



(14) 



where 5 is a three-vector consisting of the usual three Pauli matrices. 

A particle's helicity is defined as the projection of its spin on the direction of its motion. Since helicity 
of a free particle is conserved, we can choose to describe a scattering process in the basis of one-particle in- 
and out-states of definite helicity. For electrons and positrons, these states are ), where the superscript 
denotes the particle's electric charge and the subscript h = ±1/2 its helicity: 

p-S 



(15) 



States of positive helicity |e+) are often referred to as "right-handed", while states of negative helicity \et) 



%±>. 
ight-hai 

are called "left-handed". We will use this nomenclature, and replace the subscripts H — > r, > I. 

In the limit of zero electron mass, the Weyl fields cl and e R are completely decoupled in the La- 
grangian (13). After quantization, each field contains creation/annihilation operators for states of specific 
helicity only: for example, 



|0> 



l c i / : 



£l\0) ~\e+), 
e« |0> ~ |e+> , e R \0) 

Note the relation between helicity of a state and chirality of the field creating it: 

Particle : helicity = chirality ; 
Antiparticle : helicity = — chirality . 

It follows that helicity eigenstates are simply the solutions to Weyl equations of motion: 

i^d^ R = 0, i^d^ L = . 

Explicitly, the solutions have the form (up to normalization constants) 



-ip-x 



e + „ e +ip - x t 



where 



exp 



—a ■ uj 
2 



il = exp 



(16) 



(17) 



(18) 



(19) 



(20) 



where uj is the rotation from the +z direction to the direction of the momentum p. It is straightforward to 
obtain the Feynman rules in two-component language: For example, an incoming electron line of helicity h 
gives V 2E£h , an incoming positron of helicity h gives \/2E£_ h , the e~e~^j vertex is iea^, the e,~~e+7 vertex 
is iea^, etc. Note that e~e+ and eje^ vertices do not exist, since there is no coupling between and e R 
fields in the Lagrangian. 

The same construction describes the electromagnetic interactions of any other fcrmion, as long as its mass 



can be neglected: in particular, it can be applied to the e 



h+ij,- 



scattering in the limit \fs m M . 



Using the Feynman rules above, and the standard photon propagator, yields the helicity amplitudes 



M(e-e+ 
M(e^e+ 



A* ; Mr ) 

n~ni) 



M(e; 
M(e~ 



Hrl4) 



-e 2 (l-cos<9). 
-e 2 (1 + cosfl) . 



(21) 
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with all other helicity configurations giving vanishing amplitudes. Using Eq. (10), this yields the differential 
cross section 



da 



ira 



1 + cos^ 



dcos9 2s 

where we introduced the fine-structure constant a — e 2 /(4ir). 



(22) 



2.1.2. Scalar Muon 

It is instructive to repeat the above calculation replacing the muon with a scalar (spin-0) particle jl of the 
same mass and electric charge. (Such "scalar muons" , or smuons, are actually predicted by supersymmetric 
theories.) The helicity amplitudes (in the limit y/s 3> mp) are 

M(e^e+ -> fi~fl + ) = M(e~el -> = -e 2 sin0. (23) 

Homework Problem 1: Derive the amplitudes (23). 




Fig. 2. Angular distributions for muons (solid line) and scalar muons (dashed line) produced in e + e collisions. (Normalization 
is arbitrary.) 



This yields the differential cross section 

= — sin 2 6 . (24) 

ricos6> 4s v ; 

The angular distributions predicted for muons and smuons are plotted in Fig. 2. Comparing these distribu- 
tions with data allows one to distinguish between the two possible spin assignments, providing a determi- 
nation of the muon spin. (Note that spins of elementary particles are not directly observed by the detector, 
so spins can in fact only be inferred by indirect means such as this.) No evidence for spin-0 muon has been 
seen in the data, indicating that supersymmetry, if it is indeed a symmetry of nature, must be broken to lift 
the muon-smuon mass degeneracy. 

The most striking feature of the distributions in Fig. 2 is that the smuons cannot be produced in forward 
and backward directions (cos 8 = ±1), while muons can. This can be easily traced to conservation of angular 
momentum. The z-component of the total angular momentum of the colliding system is J z — ±1; all helicity 
amplitudes corresponding to J z — vanish (see Fig. 3). Since J z is conserved, the same must be true in the 
final state; but for scalar particles, J z can only be contributed by the orbital angular momentum L z . When 
particles move along the z axis, L z — 0, and angular momentum cannot be conserved; thus, scattering along 
this direction is forbidden for smuons. 



August 26, 2010 0:6 



WSPC - Proceedings Trim Size: 9in x 6in tasi09'ws'v2 



,/ 2 = -l 



Jz = +1 



Fig. 3. Helicity configurations of the e+e pair leading to non-vanishing tree-level scattering amplitudes for muon and smuon 
production. Thin arrows show the direction of particles' momenta, while thick arrows indicate their helicitics. 



2.1.3. Including the Z Exchanges 

Apart from the somewhat more complicated form of matrix elements, two conceptually important new 
features arise when the diagram involving the Z boson in Fig. 1 is included. 



nV) 




Fig. 4. Vacuum polarization corrections to the Z propagator. Each circle denotes the sum of one-loop diagrams with quarks, 
lcptons, W and Higgs bosons circulating in the loop. 



First, the Z propagator is proportional to (s — Mf) -1 , which is infinite when ^fs = Mz (as was the 
case at LEP-1 and SLC). This infinity is in fact a feature of leading-order ("tree-level") perturbation theory, 
and is automatically removed when radiative corrections ( "quantum loops" ) are included. More specifically, 
consider the vacuum polarization diagrams shown in Fig. 4. They can be included by replacing 

i i 

(25) 



P 2 - M\ p 2 - M| - n(p2) 

in the Z propagator, where Tl(p 2 ) is the one-loop vacuum polarization. We introduced the notation Mz,o 
for the bare (Lagrangian) Z mass, to distinguish it from the physical (pole) mass Mz- Note that to derive 
Eq. (25), an infinite series of diagrams in Fig. 4 had to be included and summed; this is not surprising, 
since the infinity in the tree-level calculation should be interpreted as a sign of a breakdown in perturbation 
theory. The physical Z mass is the solution to Mf - Mf - Re II (Mf) = 0. For p 2 w Mf, we can expand 

p 2 - M 2 Zfl - n(p 2 ) w p 2 - M| |0 - ( Re n(Mf) + ^Re I1(M 2 ) (p 2 - M 2 Z )\ + i Im n(M 2 ) 

1 + A Re n(M|) ) (p 2 - M%) + i Im U(M 2 ) 

= Z-\p 2 - Mf) + % Im n(M|) , (26) 

where in the last line we introduced the field strength renormalization factor Z. Note that the propagator 
no longer blows up at s — Mf , as long as the imaginary part of the vacuum polarization is non-zero. This is 
in fact guaranteed by the optical theorem, which relates Im IT to the total decay width T z of the Z boson: 

Imn(Mj) = ~ £ / ^±.^L±. m .Z ^ HM' . -^M.T„ (27) 

/l>/2 
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where the sum runs over all possible two-body decay channels of the Z. To summarize, the Z pro pagato r 
near the pole (in the 't Hooft-Feynman gauge) can be approximated by 

P 2 _ Mf + iMzTz ■ [ 0) 

If the diagram with the photon is ignored (which is a reasonable approximation near the Z pole), the cross 
section is proportional to 

2 1 

(29) 



1 



-Mf 



iM z T z 



Ml? 



M\T\ 



(s - ±r±z) 

This is the famous Breit- Wigner function. Since the above discussion did not depend on any specific features 
of the Z boson, the same function can be used to approximate the dependence of any cross section on s in 
the neighborhood of a resonance of a given mass and width. This behavior is in excellent agreement with 
experiment, as shown in Fig. 5. (The data is in fact precise enough that the contribution of the photon 
exchange diagram must be included even at the Z peak; this has been done in the figure.) 




90 92 
Vs [GeV] 




90 92 
Vs [GeV] 



Fig. 5. Cross section (left panel) and forward-backward asymmetry (right panel) of muon-pair production at LEP, measured 
by the L3 collaboration, compared with the SM prediction (solid lines). From Ref. 2 



The second interesting feature brought in by the Z is violation of two discrete symmetries, charge con- 
jugation (C) and parity (P). In the two-component language, Z coupling to electrons has the form 



with 



9l = 



-77 + < 



9R 



(30) 
(31) 
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where s Wl c w are the sine and cosine of the Weinberg angle (s^ = 0.231). The fact that gL ^ ga signals 
C and P non-conservation. The Z coupling to the muon has an identical structure, and the Z-exchange 
contribution to the helicity amplitudes for e + e~ — > n + is given by 

M{e~e+ -+ = + cos6») f B w(s) , 

M{e~ef -> n'lif) = -g 2 R {l + cos6>) f B w(s) , 

A4(e;"e+ ->• = M{e~ef -> //j>+) = -9l9r{1 - cosO) f B w(s) , (32) 

where 

fBw{3) = s - Ml + iM z Y z ■ (33) 

Summing and averaging over spins (and ignoring, for simplicity, the photon exchange contribution) yields 
the cross section 

^ = J|gf M+9 4 n) (1 + costf) 2 + 2.^(1 -cose?) 2 ] . (34) 

The cross section contains a term proportional to cos 0, so that the angular distribution is no longer symmetric 
under the parity transformation 9 — > it — 9. A useful way to quantify this is to introduce the forward-backward 
asymmetry via 

cr(cos# > 0) - cr(cos6> < 0) 
FB = cr(cos 9 > 0)+a(cos9 < 0) ' ( ^ 

It is easy to show using Eq. (34) that Af B oc (g\ — g^) 2 . Once the photon exchange diagrams are in- 
cluded, A FB becomes a function of s, since the relative contributions of the P-conserving photon exchange 
and P-violating Z exchange depend on s. Experimental measurement of Afb at LEP, along with the SM 
expectation, are shown in Fig. 5. 

2.2. Initial-State Radiation 

Not all collisions occurring at an electron-positron collider with the nominal c.o.m. energy ^/s in fact have 
that energy; in some cases, the energy is lowered by an emission of a photon (or multiple photons) by the 
electron, the positron, or both, just ahead of the collision. This process is called initial state radiation, and 
it is worth a closer look: In addition to being important in its own right, it serves as a portal into the topic 
of treatment and interpretation of infrared (soft and collinear) divergences in applications of quantum field 
theory to collider physics. These divergences are behind many of the most challenging issues in the field, 
from both conceptual and technical standpoints. 

Consider a simple extension of the muon pair-production process: the reaction in Fig. 6, e+e~ — > fi + fi~-f. 
We will first focus on initial-state radiation (ISR), described by diagrams (a) and (b); we will consider 
diagrams (c) and (d), describing the final-state radiation (FSR), in Section 2.3. Naively, this process should 
be much more rare than the muon pair-production: 

02^3 ~ —CT2^2 ~ 1^0-2^2 ■ (36) 

However, this estimate is missing an important point. Consider, for example, the diagram 6 (a). The matrix 
element is proportional to the electron propagator: 

""' h-tf-; -^; 1 <37) 

Neglecting the electron mass, the electron and photon four-momenta can be written as 

PA = (E,0,0,E); P T = {zE,p ± ,iJz*E*-v i ± ), (38) 
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e+ 





(c) (d) 



Fig. 6. Leading-order (tree-level) Feynman diagrams contributing to the process e + e — > fi + fi ■y. Diagrams (a) and (b) 
describe initial-state radiation (ISR), while (c) and (d) describe final-state radiation (FSR). 

where we introduced the photon energy fraction z = E^/E and the transverse momentum p±. Then, 



There are two limits in which the dot product vanishes, so that the matrix element in Eq. (37) blows up: 



There is also a double-singular region where both conditions are satisfied. 

What is the meaning of these singularities? First, we should note that both are artifacts of approximations 
we have made: the collinear singularity is in fact removed by including the mass of the electron, while the soft 
divergence in the tree-level diagram cancels an infinity encountered in the one-loop correction to the QED 
vertex function. So, the full theory gives finite answers. Still, the apparent divergence does have physical 
meaning: it signals a large enhancement in the probability of emission of soft and collinear photons, above 
the naive perturbative expectation of Eq. (36). This enhancement is even more pronounced in the case of 
gluon emission in QCD, where the coupling constant is larger. Understanding this effect quantitatively is 
crucial for making successful predictions. 

In an experiment, very soft and very collinear ISR photons are unobservable. Photons emerging at a small 
angle to the beam line do not get detected, since even the most hermetic detectors must have small holes 
around the beam line to let the beams in. Soft photons leave very small energy deposits in the calorimeter, 
which are drowned by noise, e.g. from thermal effects in readout electronics. Thus, experimental conditions 
set minimal values of energy and angle for which the photon is registered. To simplify the discussion, let 
us assume that a photon must have a minimal transverse momentum Q to get detected. In other words, 
if p± > Q, we register the event as a 2 -> 3 reaction, whereas if px < Q, the photon is not observed and 
the event gets recorded as an ordinary 2 — > 2 muon pair-production. The "observable" 2 — » 3 cross section 
defined in this way is automatically finite; the soft and collinear singularities are hidden in the correction to 
the 2 — > 2 cross section. Let us estimate this correction, focusing for concreteness on the collinear singularity. 

The matrix element corresponding to the diagram 6 (a) has the form 




(39) 



Pa ■ P-y — > when z — >• "soft singularity" , 

or ^— — > "collinear singularity" . 



(40) 




2p A ■ p y 



(41) 
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where the dots in square brackets denote "the rest of the diagram" (in this case, the photon propagaloi and 
the muon Dirac string; however, as we're about to show, the ISR calculation is independent of the particular 
structure of those terms). In the small-pj_ limit, p 7 w zpa up to terms linear in p±. The momentum flowing 
through the electron propagator is then pa —p-y ~ (1 — z)pa, so that (pa —p-y) 2 ~ - the electron propagator 
is nearly on-shell! This suggests replacing the numerator of the propagator with a spin sum: 

^-^E" S ((1 - z) PA ) u s ((1 - z)p A ) • (42) 

s 

The matrix element becomes 

ka - u^s uS (C 1 _ z )Pa)u s ((1 - z)p A ) , Q , a , . r , 
M 2 ^ 3 « v B Y— - 2 (p2 /2z) ' (e7 } WA £ " (P7) X [ ' ' - ] 

= -4 -EH^^-^^^M ' [wV((1-^a)] x [...] 

= -4 •E^-%e- 7 W->'^2((l-^A,P B ^Pl,P 2 ) , (43) 

where we defined the splitting amplitude 



C%e- 7 W = eu'(p A (l z))ru A e a M . (44) 

The 2 — > 3 matrix element is a product of the splitting amplitude, which describes collinear photon emission, 
and the 2 — ► 2 matrix element, describing muon pair-production in the electron-positron collision, with the 
electron energy reduced by photon emission. This behavior of the matrix element is called factorization. It 
is clear from the derivation that the splitting amplitude is universal: it does not depend on what particular 
reaction the electron enters after the photon has been emitted. The derivation can also be easily generalized 
to include soft photons. Factorization of matrix elements in soft and collinear limits is a general property of 
QED. It holds even in processes where a photon can be emitted off virtual particles in addition to external 
legs, since the diagrams with photon radiation off internal propagators are non-singular in the collinear and 
soft regimes, and their contributions are subleading. 

Physically, factorization occurs because the soft /collinear photon emission and the 2^2 scattering 
process involve different length scales. Collinear photon emission typically occurs at a distance of order m^ 1 
from the interaction point, while the 2 — » 2 scattering process involves fluctuations of fields at a much 
shorter scale of order l/y/s. The separation of scales suppresses quantum interference between the two, so 
they can in fact be considered as independent, sequential events, and total probability (proportional to the 
cross section) is simply a product of two probabilities. We will see other examples of such factorization later 
in these lectures. 

Returning to our derivation, it can be easily shown that collinear photon emission does not change the 
electron's helicity. The splitting amplitudes are only non-vanishing when s = h, where h is the helicity of 
the incoming electron: 

A4 2 -3 W "4 ■ S e-\e-,( Z ) ■ ^ 2 ((1 - Z) PA ,PB Pl,P 2 ) ■ (45) 
P± 

Plugging this form into the master formula for the differential cross section, Eq. (7), yields 

x cfflicffla 17^2^2 ((! - z)p A ,PB -> Pi,P2)\ 2 ■ <5 (4) {p A + Pb - P 1 -Pi -Pi) , (46) 



where we defined 



« - w?w <47) 
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2e 2 p\ 
z(l-z) 



l + (l-z) 2 



(48) 



independent of h. 



Homework Problem 2: Derive the formula (48). 



Using this result, Eq. (46) can be rewritten as follows: 

dcj-2 , «dn 7 ( ~— ) |i + (i-z) 2 ] 



^— — ^dilicffla \M^ 2 ((1 - z) Pa ,pb -S- Pi,p 2 )\ 2 5 (4) ((1 - z)p A +p B -p! -p 2 ) . 



2s(l - z) 



(49) 



The second line of this equation is precisely the differential cross section for 2 — > 2 scattering (in our example, 



) at a reduced center-of-mass energy — z)s. Thus, we can simply write 

(1-z) 2 ] .rfa 2 L ° 2 ((l-z) S ), 



do-' 



2->3 



fit 1 



(50) 



where the superscript "LO" serves as a reminder that this cross section is computed at tree level, or leading 
order in perturbation theory. Factorization, which was previously shown at the level of matrix elements, 
persists at the cross section level: The cross section is a product of the cross section for the reaction with no 
ISR photon, but at an energy reduced by the photon emission, and a universal factor describing the photon 
emission. Let us now work out this factor. In the collinear limit, the photon phase space element is 

d 3 p 7 1 d 2 p± ■ d(zE A ) 1 p±dp± dz 

(2tt) 3 2~K t = (2tt) 3 2zE a = 8tt 2 T ' 



rfn 7 = 



so that 



da 2 ^3 w — 



a dp± 



dz 



l + (l-z) 2 



(51) 



• da\% ((1 



z)s). 



(52) 



2ir p± 

To compute the correction to the observed 2^2 reaction rate, we need to integrate over the part of the 
phase space where the photon is unobservable, namely, p± € [0, Q}. This integral is logarithmically divergent 
at the low end. As already mentioned, this divergence in QED is not physical, but is due to the fact that 
we've set the electron mass to zero throughout the calculation. Restoring this mass is equivalent d to imposing 
a lower cutoff of m e on the p± integral in (52), which regulates the divergence: 



dz 



a 
27 



l + (l-z) 2 



log^- • do\% ((1 - Z) 8 ) . 



(53) 



This formula has a simple physical interpretation. A supposedly monochromatic beam of electrons of energy 
E = y/s/2 in fact contains electrons that lost some of their energy due to ISR, but were not deflected since 
the emitted photon was collinear, and so may collide with the positrons and initiate reactions such as muon 
pair-production. The probability of finding an electron carrying energy between xE and (x + dx)E, where 
x € [0, 1), can be read off from Eq. (53): since x = 1 — z, we obtain 

"I + x 2 ~ 



f e (x)dx = — 



l-x 



log — dx . 



(54) 



d The equivalence works in the limit Q S> m e ; that is, the formula (53) captures the leading (logarithmic) behavior in this limit, 
but not the subleading corrections. 
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Of course, the total probability of emitting an ISR photon is less than one, so f e (x) must also cunLain 
a term proportional to 5(1 — x). In addition, since the collinear ISR photons travel in essentially the same 
direction as the original electron, they may also collide with the positrons (initiating, for example, a Compton 
scattering process 7 + e + — > 7 + e + ), so it is quite reasonable to consider them a part of the beam as well. 
The probability of finding a photon with energy between zE and (z + dz)E is given by 

f 7 (z) = / e (l - z). (55) 

Since both electron and positron beam contain photons, there is even a possibility of a photon-photon 
collision: Fig. 7 shows the cross section for 77 — > measured by the L3 collaboration at LEP. 




(GeV) 

Fig. 7. Comparison of data from the L3 experiment at LEP-2 (points) and theoretical prediction (solid line) for the cross 
section of 77 — > as a function of 77 center-of-mass energy. From Ref. 3 

The functions f e and / 7 describe the composition of the electron, in a way similar to how the parton 
distribution functions (pdf's) of the parton model describe the composition of the nucleon (more on pdf's 
in Lecture 3). When ISR is included, an "electron beam" in fact consists of electrons and photons. The 
probabilities of finding electrons and photons in the beam depend on Q: This should not be surprising, since 
Q sets the boundary in pt below which the ISR photons would not be observed and so are included in the 
beam. In particular, if all photons were observed, corresponding to Q = m e , we would get f e = 6(1 — x), 
f-y = 0, consistent with the naive picture of the beam. 

The distribution f e (x, Q) is singular when x — > 1, while / 7 (z, Q) blows up when z — > 0. This is the region 
where the collinear ISR photons are also soft, i.e. overlapping soft and collinear divergence. As I already 
mentioned, this divergence is canceled when next-to-leading order (NLO) contribution to the 2 — > 2 cross 
section, containing the interference term between the tree-level and one-loop matrix elements, is taken into 
account. A convenient way to see this cancellation is to introduce a fictitious photon mass /i, which regulates 
the z — » divergence in (52) in the same way as the electron mass regulates the p± — »■ singularity. The 
NLO contribution to the 2^2 cross section also has logarithmic dependence on /z, coming from an infrared 
singularity in the one-loop correction to the matrix element which is regulated by the photon mass. The 
dependence on /j, cancels in the sum of ISR and NLO cross sections, giving a finite and well-defined answer 
for f e and / 7 . However, the probability of emitting a single soft photon is typically not small, indicating that 
processes with emission of multiple soft photons are important. The reason for this is that while formally 
each extra photon in the final state costs a factor of a/n ~ 1/300 in the cross section, in reality collinear 
photon emission is only suppressed by (a/n) log(Q/m e ), with an additional logarithmic enhancement in the 



c This is another example of a breakdown in perturbation theory in a particular region in phase space, similar to its breakdown 
in the neighborhood of an s-channel resonance discussed above. 
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soft-collinear double-singular region. The logarithmic factors can be large, slowing down the converg ence o f 
the perturbation theory. To obtain accurate predictions, terms of higher order in a must be included, at 
least those containing the highest power of the large logarithms possible at each order. Luckily, this problem 
turns out to be tractable. To all orders in a, leading-logarithm accuracy, the distribution functions obey a 
system of integro-differential equations known as Gribov-Lipatov (GL) equations. For example, 

^Sf = ?f f [W-)(/.-(f0)+^(f,9)) + i^WA(f.0)]. m 

where f e + is the positron distribution function. (Positrons must be present in the "electron beam" due 
to processes such as — > e~e~e + , which appear beyond the leading order in perturbation theory.) The 
splitting functions are given by 

_ l + (l-z) 2 

Z 

P^ = ~S{\-z). (57) 

Note that P e ^ 7 is essentially the leading-order ISR photon emission probability computed above, while 
P 7 ^ 7 is proportional to a delta function, reflecting the impossibility of a photon splitting into a photon pair 
at leading order in QED. (The photon can, however, split into an e+e~ pair, hence the coefficient of the 
delta function is not unity.) The physical meaning of this equation is as follows: if we shift Q — > Q + SQ, 
the photons with p± G [Q, Q + SQ] that were previously counted as "detectable" should now be counted as 
"part of the beam" and included in / 7 . The beam at Q consists of electrons, positrons, and photons; each 
of these particles, if its energy is above xE, can emit an ISR photon with energy xE and p± £ [Q,Q + SQ]. 
The right-hand side of the equation simply sums up the probabilities of such emissions. Solving Eq. (56), 
together with the corresponding equations for f e ± , with the boundary conditions 

f e -(x,Q = m e ) = S(l - x),f 7 (Q = m e ) = f e +(Q = m e ) = 0, (58) 

gives the desired distribution functions. For more details, as well as for a derivation of the GL equations, the 
interested reader is encouraged to read Section 17.5 of Peskin and Schrocdcr's textbook. 1 



2.3. Final State Radiation 

Final-state radiation diagrams, Fig. 6 (c) and (d), can be treated in exactly the same way as the ISR. Again, 
the cross section is dominated by photons that are either soft or collinear with the muon or antimuon. One 
slight difference experimentally is that collinear (but not soft) FSR photons can be observed, so instead 
of using p± to separate 2 — > 3 and 2^2 events, we simply assume that the photon gets detected once 
its energy is above i? m i n . Following the same steps as in the derivation of Section 2.2, we obtain the total 
observable 2^3 cross section: 

a^+AO) = - log log • <t L0 (// + /O , (59) 

where we assumed y/s 3> m M and y/s 3> -Emm, and ignored terms not enhanced by large logarithms. The 
factors of log -gf— and log come from the soft and collinear singularities, respectively; their product is 

min M 

referred to as the Sudakov double logarithm. To obtain the correction to the observed 2^2 rate, we must 
combine the contribution from FSR with E 1 < E m - ln with the NLO correction to the 2—^2 cross section. 
The result (again in the leading-log approximation) is 

a° hs (n+»-) = (l - - log -|- log JL) • <7 LO ( M + M -) . (60) 
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Notice that, up to terms with no large logarithms, cr obs (/i + ^~) + o-(n + fi~j) — er LO (/i + /i~). This result is 
independent of E min . This observation will play an important role when FSR in QCD (gluon radiation off 
quarks) is considered below. 

2.4. Hadronic Final States 




Fig. 8. Leading-order (tree-level) Feynman diagrams contributing to the process e + e — > qq. 

Since quarks are electrically charged and couple to the Z, quark-antiquark pairs can be produced in e + e~ 
collisions in the same way as muon pairs, see Fig. 8. A qualitatively new feature in this case is that the quarks 
themselves are not observed: they turn into color-neutral hadrons, which are then detected by the tracker 
(if electrically charged) and the hadronic calorimeter. A complete theoretical description of the scattering 
process must include a model of quark hadronization. Since hadronization is inherently non-perturbative, 
involving momentum transfers comparable to the QCD confinement scale Aqcd ~ 200 MeV, it is very 
challenging to analyze this part of the reaction from first principles, and in fact all existing approaches rely 
on phenomenological models. It is hard to quantify their accuracy. In general, the most useful theoretical 
predictions involving hadrons are those that are least dependent on the details of hadronization. This typically 
involves sums or averages over a large class of final states. Observables involving such sums are often referred 
to as inclusive, while those where a particular final state is specified are called exclusive. An extreme example 
of a fully inclusive observable is the total e + e~ — > hadrons rate considered below: At collision energies large 
compared to Aqcd (which we will assume throughout this section), it can be predicted, up to corrections 
suppressed by powers of Aqcd / v^i with no recourse to a hadronization model. On the other extreme are 
fully exclusive quantities, e.g. cr(e + e~ — > 2tt + 2K + 3p), which can only be calculated within a hardonization 
model. Trustworthy quantitative predictions for such observables are currently out of reach. Luckily, there 
exists a large class of useful observables that are sufficiently inclusive to allow for precise predictions, and at 
the same time carry much more information than the total hadronic rate. We will consider an example in 
Section 2.4.3. 

2.4.1. e + e _ — > Hadrons 

The simplest inclusive observable is the total hadronic cross section, a(e + e~ — > hadrons). If collision energy is 
well above Aqcd , it is reasonable to assume that there is no quantum interference between the short-distance 
process of quark pair creation (characteristic length scale 1/^/s), and the long-distance hadronization process 
(typical scale ~ Aq CD ) f ■ ^ n °th er words, 

a(e + e~ —> hadrons) = cr(e + e~ — > qq) x Prob (qq —> hadrons) , (61) 



This is another example of factorization due to separation of scales, see comments under Eq. (44). 
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up to corrections of order Aqcd/v^ an d ct s - Since each quark must hadronize in one way or another, I he 
probability factor is equal to unity. A simple calculation then gives (for collision energies well below mz) 

a(e + e~ -> hadrons) = 3 • Q q ■ a(e + e~ -> m + aO , (62) 

9 

where Q q are the quarks' electric charges, and only quarks with m q -C should be included in the sum. 
(In threshold regions, a/s f=a 2m q , the cross section behavior is somewhat more complicated due to quark 
mass effects and the presence of bound states.) It is customary to present experimental data on this process 
in terms of a dimensionless ratio 

a(e + e~ — > hadrons) 
R = -i_ J. . 63 

cr(e+e~ -> fi + fi~) 

At energies well below the Z peak, and away from thresholds, the formula (62) is in excellent agreement 
with experiment, see Fig. 9. This gives firm support for the parton model and the factorization assumption. 
The calculation can be easily extended to include Z exchanges. This result is in good agreement with data 
up to the highest energies reached at LEP-2, see Fig. 10. Radiative corrections must be included to match 
the high experimental precision available for y/s « Mz (see section 2.4.3). 
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Fig. 9. Total hadronic e + e scattering cross section, normalized to rj(e + e — > fi + fJ. ), as a function of center-of-mass energy, 
at energies well below the Z peak. The solid blue line is the SM prediction. From Ref. 4 



2.4.2. Hadronization 

To predict cross sections for specific hadronic final states (e.g. 7r + 7r~ or K + K~tt tt ), one must model 
the hadronization process. In practice, this is done by using a Monte Carlo (MC) generator, a computer 
program implementing one of the available phenomenological hadronization models. Popular hadronization 
models include the string fragmentation model, implemented in PYTHIA, 6 and the cluster hadronization model, 
implemented in HERWIG. 7 The models are quite complicated, and contain a large number of free parameters 
that must be fitted to data. A detailed description of hadronization is beyond the scope of these lectures. 




Still, just to give a flavor of the physics involved, let me briefly consider the basic picture behind the string 
fragmentation model. 




Fig. 11. String fragmentation model of hadronization. (a): Quark-antiquark pair, connected by a color flux tube, immediately 
after the collision, (b): After the first string breakdown. Filled circles denote primary partons, while open circles denote secondary 
partons. Primary parton momenta are i/s, while the typical secondary parton momenta are of order Aqcd <C yfs- 



Consider a quark-antiquark pair created in an e + e collision. Let us refer to these q and q as "primary 
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partons" . To hadronize, the primary quark must form a bound state with an antiquark; however, the p rimar y 
partons are rapidly moving away from each other. Their invariant mass yfs is large compared to the binding 
energy <~ Aqcd, so they typically do not bind. Instead, quark hadronization proceeds by a spontaneous 
creation of new "secondary" quark- antiquark pair(s) from the vacuum, and formation of bound states between 
the primary quark and a secondary anti-quark. The spontaneous pair-creation process is possible because the 
primary partons create color (gluon) field in the surrounding space; it is this field that supplies the needed 
energy. Color field lines originate at the quark and end on the antiquark. In the string model, the color field 
immediately after the collision is modeled as a color flux tube, or string, connecting the primary partons: see 
Fig. 11. The string tension is of order Aq CD ~ 1 GeV/fm, and its width is of order Aq CD , corresponding 
to typical transverse momenta <~ Aqcd- As the primary partons move apart after the collision, the string 
between them is stretched, and at some point it breaks down. A secondary qq pair is spontaneously created 
at the breakdown point, so that the original string splits into two string fragments. The secondary pair has 
typical momentum (both along and transverse to the string) of order Aqcd, so that the invariant mass of 
each string fragment after the breakdown is of order y/E q h.QCT> ~ s^Aqcd ^ Each string fragment 
then breaks down again, further reducing the invariant mass. The process is iterated until the invariant 
mass of all fragments is about 1 GeV. At this point, each string fragment is associated with a meson. The 
probabilities for a string to turn into a particular meson are phenomenological parameters, to be fitted to 
data. A similar, though a bit more complicated, picture is used to describe baryon formation. 

It is clear that the above model is heuristic, and a large number of input parameters is needed: string 
tension, breakdown probability, secondary pair momentum distributions, relative probabilities of nucleating 
uu, dd and ss pairs, etc. A reasonable agreement with data can be reached by adjusting these parameters; 
however, it is difficult to quantify how uncertain the model predictions really are, especially for highly- 
exclusive observables. 

2.4.3. e+e- -t Jets 

Since secondary quarks and antiquarks have typical momenta of order Aqcd <C y/s, most of the momentum 
carried by the hadrons is due to the primary partons. This can be seen explicitly in the string model 
outlined above, but is in fact independent of the details of the hadronization model. The hadrons come 
out in two collimated jets, one approximately collinear with the quark and the other with the antiquark. 
(Experimentally, a "jet-finding algorithm" is used to find clusters of approximately collinear hadrons in 
the event and identify them as jets, and the precise definition of a jet depends somewhat on the details 
of this algorithm.) The sum of the four-momenta of all hadrons in each jet is equal to the momentum of 
the corresponding primary parton. Corrections to this expectation, and the opening angle of each jet, arc 
of order 8 Aqcd/\/s- While these corrections depend on the hadronization model, they are small. Ignoring 
them, jet production cross sections and kinematic distributions are identical to those of the primary partons, 
which can be calculated within perturbation theory. These cross sections and distributions provide a wide 
class of observables that can be used to test the SM. 

As an example, consider the angular distributions of jets produced in e + e~ — > qq. The differential cross 
section da/d cos 9 can be easily computed: at tree level, it has the same form as the muon cross section in 
Eq. (34), but with 




(64) 



g This estimate applies only to effects of hadronization. As we will see below, there are additional corrections from higher-order 
perturbative QCD; however, these are still small, and in addition are (at least in principle) fully calculable. 
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where the superscript refers to up-type and down-type quarks. Comparing this prediction with data at th e 
Z resonance could be used to test the SM prediction for quark couplings to the Z boson, in the same way 
as the measurement of the muon forward-backward asymmetry discussed in Section 2.1 tests the muon cou- 
plings. However, there is an additional complication: Jets originated by u, d, s quarks and their antiparticles 
are indistinguishable experimentally, since it is impossible to determine which of the partons making up 
the hadrons in a given jet was the primary. This ambiguity completely washes out the forward-backward 
asymmetry, which has opposite sign for particles and antiparticles. The situation is better if the primary is 
a heavy quark, c or b: since the probability to nucleate heavy quark pairs from the vacuum is very small, 
the jet typically contains just a single heavy quark. This heavy quark can be detected using the fact that it 
travels a small but macroscopic distance before decaying (cr w 500 /jm for b, 100 — 300 /im for c), producing 
a secondary vertex in the inner tracker detector. If the decay is semileptonic, the charged lepton inside the 
jet can serve as an additional signature. A combination of these techniques is known as flavor-tagging: For 
example, a jet with an identified b quark is said to be "b-tagged" . In addition, the charge of the lepton from a 
semileptonic decay can be used to infer the charge of the heavy quark (e.g. b — > l~ vs. b — > £ + ; make sure you 
understand why!), giving a "charge tag". This allowed for a direct measurement of the forward-backward 
asymmetry in e+e" -> bb at LEP-2 and SLC. The result is A^ B = 0.0992 ± 0.0016, about 2.5 standard 
deviation below the SM prediction. This is the largest deviation from the SM among precision electroweak 
observables, although not large enough to be taken seriously as an indication of new physics. 




Fig. 12. Feynman diagrams contributing to the process e + e — > qqg at the tree level. 



At leading order in the QCD coupling constant a s , all hadronic events in e + e~ collisions result in two 
jets. At the next-to-leading order, an additional gluon can be emitted, see Fig. 12. While formally suppressed 
by a s /7r ~ 1/30, this process is enhanced by singularities for soft and collinear gluon emission, analogous 
to those encountered in the photon FSR analysis of Sec. 2.3. The emitted gluon hadronizes, creating an 
additional jet. The experimental signature of this jet depends on its direction. If the gluon is approximately 
collinear with the (anti-)quark, the jet-finding algorithm will merge the hadrons originating from the (anti- 
)quark and the gluon into a single jet, identifying the final state as 2-jet. If the gluon is emitted at a large 
angle, the extra jet will be reconstructed separately, resulting in a 3-jet event. Thus the 2-jet and 3-jet rates 
depend strongly on the details of the jet-finding algorithm, and the value of parameters used to define a jet. 
However, just as in the case of photon FSR, the sum of the two rates is completely independent of these 
details. A complete NLO-QCD calculation gives h 

^(e+e" ->> 2 jets) = a LO (e+e- -> qq) ' 



1 

47T 



(65) 



Comparing this result with the cross section at the Z resonance precisely measured by LEP experiments 
yields one of the most precise determinations of the QCD coupling constant: a s (Mz) = 0.123 ± 0.004. 



h A reader interested in the derivation of this formula is encouraged to work through the Final Project at the end of Part I of 
Pcskin and Schroeder's textbook. 1 
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Numerically large value of a s and collinear enhancement mean that most jets contain multiple c ollinea r 
gluons in addition to the primary parton. These gluons can in turn split into qq or gg pairs, which can emit 
additional gluons etc. This phenomenon is known as parton showering. The internal structure of the jets, i.e. 
observables such as the jet opening angle, is dominantly determined by the parton showering, which plays a 
more important role than hadronization. In addition, depending on the parameters used to define jets, parton 
showers may contribute significantly to the observed n-jet rates for n > 2. Parton showering can in principle 
be completely described in terms of perturbative QCD, although resummation of large logarithms to all 
orders in perturbation theory is necessary. In practice, most phenomenological and experimental analyses 
rely on the numerical implementation of parton showers in Monte Carlo programs such as PYTHIA and 
HERWIG. 

3. Hadron Collisions 



Pi / \ XiPi \J x 2 P 2 



Fig. 13. Parton model picture of a hadron-initiated process p(P\) + p(ft) — ► Y + X. 



Theoretical analysis of hadron collisions (pp or pp) proceeds along the same lines as for e + e~. The 
main new difficulty is that the initial state consists of composite particles. According to the parton model, 
production of final states with total invariant mass large compared to Aqcd is initiated by a pair of partons, 
with the rest of colliding hadrons serving as spectators, as represented graphically in Fig. 13. Quantitatively, 
for a proton-proton collision, this picture implies that 



da{p{P 1 )+p{P 2 )^Y + X) = f dx x [ dx 2 fi t (xi) fi 2 {x2)da(i 1 (x 1 P 1 



)+i2{x 2 P 2 )^Y). (66) 



Here, Y denotes the high-invariant-mass final state of interest; X denotes anything else (including the re- 
hadronized remnants of the colliding protons); Pi and P 2 are the momenta of the incoming protons; x\ and 
x 2 are the fractions of those momenta carried by the partons that initiate the reaction; and the sum runs over 
the types of partons (quark and antiquark flavors and gluons). The functions fi(x) are the parton distribution 
functions (pdf's); we will discuss them in more detail below. The analogous formula for cross sections in pp 
collision is obtained by replacing fi 2 (x 2 ) — > fi 2 (x 2 ) — fi 2 (x 2 ), where g = g by definition. The parton-level 
cross section on the right-hand side, a(ii(xiPi) + i 2 (x 2 P 2 ) — > Y), is calculated within perturbation theory. 

While the center of mass of the colliding hadron pair is at rest in the lab frame, the parton frame, i.e. 
the frame in which the center of mass of the colliding partons is at rest, is moving with respect to the lab 
frame with velocity 

[3 = (67) 
Xi + x 2 

along the collision axis. (Note that the parton-level cross section in Eq. (66) can always be evaluated in the 
parton frame, since the differential cross section is invariant with respect to boosts along the collision axis.) 
Conservation of energy-momentum implies that the invariant mass of the state Y has to be equal to the 
parton center-of-mass energy, y/I, where 

S = X\X 2 S . 
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Since x\ and x-i are not the same in each event, but instead are picked from a distribution dictated by th e 
pdf 's and cross sections, the values of /3 and s in each event are a priori unknown. If the final state Y is fully 
reconstructed (all particles are detected and their energies and momenta measured), these quantities can be 
experimentally determined on an event-by-event basis. If, however, Y includes invisible particles (neutrinos 
in the SM, neutralinos in supersymmetric models, etc.), their values in each event remain unknown. This 
lack of kinematic information complicates the analysis of such final states at a hadron collider. 



3.1. Parton Distribution Functions 

Building qualitative intuition about hadron collisions requires some familiarity with the parton distribution 
functions. Conceptually, the pdf's are close cousins of the electron and photon distribution functions in an 
"electron beam", encountered in Section 2.2. Part of their job is exactly the same: They account for emission 
of collincar ISR gluons, whose transverse momentum px is too small for them to be detected individually, 
and for further splitting of these gluons into collinear qq and gg pairs. Just as in the photon ISR case, 
these effects depend on the minimal transverse momentum Q for which a gluon (or a quark from gluon 
splitting) can be resolved as an additional jet; thus, the pdf's are really functions of two arguments, x and 
Q. The dependence of the pdf's on Q is encoded in the Altarelli-Parisi (AP) equations, which extend the 
Gribov-Lipatov equations of Section 2.2 to QCD. For example, the counterpart of Eq. (56) is 

df g (x,Q) _ a s f 1 dz 
aiogQ 

where the QCD coupling constant a s should be evaluated at the scale Q. The splitting functions are given 

by 

4 rl + (l-z) 2 " 



f f [p q ^)(f,(- z ,Q)+h(- z ,Q))+P g ^)f,(- z ,Q)} , (69) 



p 

g ^ 9 ~ 3 



P'y — yy — 6 



z 



" + — +z{l-z)+( l ±- r ^\5{l-z) , (70) 



-{l-z)+ z v ' \12 18, 
where nf is the number of quark flavors with masses m q -C Q. (For derivation, as well as definition of the 
"+ prescription" used to regulate the denominator of the first term in P g ^ g at z — > 1, see Section 17.5 of 
Peskin and Schroeder. 1 ) The physical meaning is the same as for the GL equation, see comments at the end 
of Section 2.2; the only qualitative difference is the possibility of g — > gg splitting at leading order, due to 
non-Abelian nature of QCD. The key new feature is the boundary conditions at small Q. In the case of QED, 
the conditions (58) are motivated by the observation that, if all ISR photons were detected, one would simply 
end up with a monochromatic electron beam. In a hadron collision, though, even if every collinear gluon is 
detected, one still needs to specify the probabilities of finding various partons in the proton. These cannot be 
computed within perturbation theory, and at present the only viable approach is to measure them from data. 
Thus, in applying the parton model, one must first perform a fit to a number of experimental observables 
which would provide the pdf's, and then use these pdf's to predict other observables. The reason this works 
is that the pdf's are universal, that is, independent of the reaction that the partons enter. The proof of this 
universality follows along the same lines as the collinear factorization proof in Section 2.2, although details 
are considerably more involved. 

Which observables should be used to determine pdf's? They should obey three important criteria. Firstly, 
the set should be broad enough to achieve sensitivity to as many pdf's and in as broad a range of x values as 
possible (Q values are less of a concern, since evolution in Q is calculable via the AP equations.) Secondly, 
an accurate theoretical prediction at the parton level should be available. Observables used in practice are 
typically known to NLO or NNLO in the QCD coupling constant. Finally, only processes where contamination 
from new physics beyond the SM is unlikely should be used; ideally, all data included in the fit should be at 
energies well below the scale where the SM has already been extensively tested and verified. 
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Table 2. The main processes included in the MSTW-2008 fit, ordered in three 
groups: fixed-target experiments, HERA, and the Tcvatron. For each process, the 
table lists dominant partonic subprocesses, the partons whose pdf's are primarily 
probed, and the x range constrained by the data. From Ref. 8 



Process 


Subprocess 


Partons 


x range 


l± {p, n}-^f ± X 


7*5 -> g 


q,q,g 


x > 0.01 


nip -» £ ± X 


7* ri/u — » d/u 


d/u 


x > 0.01 


pp — > l^ + fJ.~ X 


mm, dd — > 7* 


q_ 


0.015 < x < 0.35 


pn/pp — ¥ u + fi~ X 


(ud)/(uu) — > 7* 


d/u 


0-015 < x < 0.35 


u(u)N -» ^-(a+)X 


W*q -> q' 


q,q 


0.01 < x < 0.5 


vN -» x 


W*s -> c 


s 


0.01 < x < 0.2 


UN -> n+fi~ X 


W*s -» c 


s 


0.01 < x < 0.2 


e ± p -y e ± X 


7*9 -> 9 


g,q,q 


0.0001 < x < 0.1 


e+p^uX 


{d,s} -> {m,c} 


d, s 


x>0.01 


e^p — > ccX 


7*c — > c, 7*<y — > cc 


c, g 


0.0001 < x < 0.01 


e ± p — » jet + X 


1*9 QQ 


g 


0.01 < x < 0.1 


pp -> jet + X 


gg,qg,qq 


g,q 


0-01 < x < 0.5 


pp-> (w ± -»■ e k v)x 


ud — > W, -> W 


m, d, m, d 


x > 0.05 


pp-> {z-> e+e-)x 


mm, — » Z 


d 


x > 0.05 



As an example, consider the recent MSTW-2008 pdf fit. 8 The reactions used in the fit, along with 
the pdf's primarily constrained by each reaction, are listed in Table 2. Most of the data comes from deep 
inelastic scattering (DIS) experiments, including neutral-current scattering ep — > eX and charged-currcnt 
reactions ep — > vX, vp — > IX. The most accurate and complete data set on ep scattering is provided by 
the experiments at the HERA ep collider. Using DIS data has an important advantage that only one of the 
reacting particles is composite. A disadvantage is that only some linear combinations of pdf's are determined. 
For example, there is no direct sensitivity to the gluon pdf, since at tree level leptons do not interact with 
gluons. Measuring quark pdf's at different Q scales provides an indirect measurement of the gluon pdf via 
AP evolution equations. This method works quite well at low x. The fit includes additional data, such as the 
Tevatron dijet rate (see Sec. 3.3), to constrain the gluon pdf at high x. Large uncertainties remain at x > 0.1, 
partly because f g is quite small in that range. Additional constraints on quark pdf's are provided by the 
measurements of electroweak gauge boson rapidity distributions at the Tevatron, as discussed in Sec. 3.2.2. 

Figure 14 shows the MSTW-2008 pdf's for two values of Q: Q 2 = 10 GcV 2 and 10 4 GcV 2 (the latter 
corresponding roughly to Q = Mz, a common choice for describing electroweak processes). Error bars are 
indicated by the width of the lines; roughly speaking, pdf's can be varied within these bands without spoiling 
consistency of the fit with the data. A few features of the pdf's are worth remembering for anyone wishing 
to build qualitative intuition about hadron collider physics. First of all, note that while in the naive quark 
model the proton is "made of u and d quarks, in reality they only dominate the pdf's at x 0.1. (The u 
and the d do carry all of the proton's electric charge, but the pdf's describe the distribution of momentum, 
not charge, among partons.) For lower x, the dominant component of the proton is the gluon, especially in 
the low-x and high-Q 2 regions. (Note that in the figure, the gluon pdf is divided by 10!) This is especially 
relevant at the LHC, where most collisions with parton center-of-mass energies up to about 1 TeV will be 
between gluons; for this reason, the LHC is sometimes called a "gluon collider". Moreover, below x ^ 0.01, 
the pdf's of all light quarks and antiquarks are roughly the same; this is because most of them come from 
gluon splittings, which are flavor-blind and produce equal number of q and q. (The strange and charm pdf's 
are suppressed due to their mass, but the effect is not strong. At Q ^> m,b, bottom pdf should be taken into 
account as well.) The low-x quarks and antiquarks are often referred to as the "sea" partons, whereas u and 
d quarks carrying x > 0.01 are called "valence quarks". Of course, the distinction between sea and valence u 
and d quarks is not precise, so these terms can only be applied in a qualitative sense. As far as Q evolution 
is concerned, for both quarks and gluons, pdf's migrate towards lower x as Q is increased. The reason is 




Fig. 14. MSTW-2008 pdf's at Q 2 = 10 GeV 2 and 10 4 GeV 2 . The width of the lines indicates the error bars. From Ref. 8 

that, as Q is increased, new partons are added that come from collinear splittings of the original partons. 
Since splitting always lowers x, this increases density at low x and suppresses it at high x. This effect can be 
clearly seen in Fig. 14. Since splitting amplitudes are proportional to the QCD coupling constant evaluated 
at the scale Q, the speed of the evolution decreases with Q due to asymptotic freedom of QCD. 

3.2. Electroweak Gauge Boson Production 

As an example of a hadron collider process, let us consider production of a single Z boson. We will first 
compute the total production cross section of the Z at the Tevatron and the LHC, and then proceed to 
discuss its kinematic distributions. 

3.2.1. Cross Section 

At leading (tree) level in perturbation theory, the Z can only be produced in qq collisions, with cross section 

<r(M^Z) = ^^^6(s-M>), (71) 

where T(Z — > qq) = I 'z ■ Br [Z — > qq) is the partial decay width of the Z in the qq channel. At the hadron 
level, this yields 

a(pp^Z + X) = ^-^f- [ dxi [ dx 2 y>/,(a;i, Q)fq(x 2 ,Q) Br (Z ^ qq) 6(x lX2 s - M%) . (72) 
A Jo Jo q 

The same formula applies to pp collisions, with the substitution 

2f q ( Xl ,Q)fg(x 2 ,Q) — > j q {x u Q)f q (x 2 ,Q) + f q ( Xl ,Q)f q {x 2 ,Q) . (73) 
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An alarming feature of Eq. (72) is that the cross section seems to depend on the scale Q. Recall th aL Q i s 
defined as the minimal transverse momentum at which an ISR parton is registered as an extra jet. However, 
we are considering the process Z + X, where X includes any hadronic activity, so events with any number of 
extra jets should be included and the cross section should be independent of Q. The reason for this apparent 
contradiction is that we used the leading-order parton-level cross section in our calculation. When radiative 
corrections are included, the parton-level cross section itself becomes Q-dependent. For example, at NLO 
in a s , the parton-level cross section contains a contribution from the process qq — > gZ, but only gluons 
with pr > Q should be included. (The gluons with pr < Q have already been accounted for by the pdf's, 
so including them again in the parton-level cross section would be double-counting.) Including the NLO 
parton-level cross section cancels the Q dependence of the pdf's up to terms of order a 2 s ; cancelling those 
terms would require an NNLO parton-level calculation, etc. The upshot is that while the true hadron-level 
cross section is Q-independent, in practice theoretical predictions always have a residual Q-dependence due 
to uncalculated higher order terms in perturbation theory. This raises two questions: How does one choose 
Q! And, since the Q dependence should presumably be considered as a systematic uncertainty on theoretical 
predictions, how does one assign a conservative but reasonable value for this uncertainty? Not surprisingly, 
there are no precise answers. Several prescriptions have been proposed based on physical arguments. As more 
and more higher-order calculations become available, these prescriptions can be tested against "data" : For 
example, if both NLO and NNLO answers are known, as is the case for pp — > Z + X, one can ask whether 
the central value and the error bar based on an NLO calculation would be compatible with the NNLO result. 
Without going into details, let me simply state the simplest prescription: For the central value, take Q to be 
the invariant mass of the final state, in our case, Mz- (This is motivated since higher-order terms involving 
collinear emission are enhanced by log(s/Q 2 ), as we saw in Section 2; the choice Q ~ Mz avoids large logs.) 
For the error bar, take the variation of the cross section as Q is varied between Mz/2 and 2Mz- This recipe 
is often used in practice, and tends to roughly agree with more sophisticated prescriptions. 

Returning to our calculation, we can easily perform one of the integrals: 

f 1 1 

/ dx 2 S(x 1 x 2 s- Ml) = Oixxs-Ml), (74) 

Jo x i s 

yielding 

a(pp^Z + X) = ^^f- f V q {xuM z )f- q {^,M z ) Br (Z -+ qq) . (75) 

3 M Z S J Ml/a x l g x i s 

Note that only partons carrying momentum above a certain threshold can participate in Z production: 

Ml _, 
Xmin = — - ~ 2.5 x 10 at the Tevatron, 
s 

w 4.0 x 10~ 5 at the LHC . (76) 



Comparing with Table 2, we confirm that pdf measurements at HERA probe the entire range of x values 
necessary for the Tevatron calculation, and almost the entire range (except for the extreme low- a; corner) 
needed for the LHC. (In practice, the cross section only receives a small contribution from x <~ £ m in, since 
this requires the other parton to be close to x = 1, where pdf's are suppressed. Most of the cross section 
is typically contributed by the region with x\ ~ X2 ~ ^min? 

although this is not a sharp peak.) The next 
step is to perform the x\ integral. This can be done numerically using the pdf's downloaded from the web 
site 9 in C++, Fortran, or Mathematica-compatible versions. To get an estimate, we can use rough log-linear 
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extrapolations of the pdf's: for example (c/. Fig. 14) 



xu(x, Mz) 



For the Tevatron parameters, this gives 



-1.2-0.91og 10 z, 
0.7, 

-0.2- 0.91og 10 x, 
0, 



dx i . f x m j r 

uixi) u 

xi v \ Xi 



10 



dx\ _, \ — ( X m {i 

u(xi) u 

xi V x i 



<x < 10~ 2 ; 
10~ 2 < x < 0.1; 
O.K x < 0.6; 
x > 0.6. 



600, 
5. 



(77) 



(78) 



So, Z production at the Tevatron is mostly due to valence quark collisions, with sea partons contributing 
only about 1% of the cross section. Plugging in the numbers gives the cross section (from uu collisions alone) 
of about 3 nb. The actual SM prediction, including all quark flavors as well as NLO-QCD corrections, is 
8.2 nb. Note that the NLO (0(a s )) corrections are not small, enhancing the cross section by about 30%. 
(This is not unique to Z production: many hadron-initiated processes receive an NLO-QCD correction of a 
similar size or even larger.) Next-to-next-to-leading order (NNLO) corrections have also been computed 10 ' 11 
and are numerically much smaller, indicating convergence of the perturbation series. The SM predictions for 
inclusive Z production, and for a very similar process pp — > W + X, are in excellent agreement with data, 
as shown in Fig. 15. 
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Fig. 15. pp — > Z + X and pp — ¥ W + X cross sections, multiplied by leptonic branching ratios, as a function of center-of-mass 
energy. Data collected by CDF and D0 collaborations at the Tevatron Run-I and Run-II, as well as earlier data from the CERN 
SPS collider, are shown. Solid lines: SM prediction. From the CDF collaboration web site. 12 



At the LHC, we obtain 



^l u ( Xl )u{^) w 6 x 10 5 , (79) 



X\ \ X\ 
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and the total cross section is about 60 nb. Z production at the LHC is dominated by sea parton col lisions , 
and the cross section is larger than that at the Tevatron, thanks to the growth of pdf 's at low x. The growth 
of cross section with \fs is due to the composite nature of hadrons, and is in sharp contrast with the decrease 
of cross sections with y/s in elementary particle collisions, see Eq. (6). 



3.2.2. Z Rapidity Distribution 




Fig. 16. Kinematics of Z — > £ + £ decay in the parton frame. 



The produced Z boson is at rest in the parton frame; in the lab frame, it has no transverse momentum, 
and is simply moving with velocity f3, see Eq. (67), along the beamline. The Z decays promptly; let us 
consider the leptonic decay Z — > , where I = e or fi. (The "Drell- Yan" process pp — » Z — » i + i ~ provides 
the cleanest signature of the Z at hadron colliders due to low backgrounds. Hadronic Z decays must be 
distinguished from a large background of 2-jet events from pure QCD processes, which is challenging.) The 
kinematics of the leptons in the parton frame is shown in Fig. 16. For each lepton, its pseudorapidity fji is 
defined by 

tank fji = ^ , (80) 
Ei 

where hats on all symbols indicate that they refer to parton-frame values of the variables. (The term "pseu- 
dorapidity" is often abbreviated to simply "rapidity" in collider physics applications, since true rapidity is 
never used; we will follow this practice below.) With this definition, 

P\\i =p±ismhfii , 

E i = \J (.hi) 2 + (P-Li) 2 = P±i COsh Vi , (81) 

where = |p_u|. Note that p±i = j5_i_2, and fji = —fj2- To find lepton momenta in the lab frame, we must 
boost by -0 Z = ff^. Defining 

r) z = tarni- 1 /?,^ 1 lo %TZJ z ' ( 82 ) 

it is easy to show that 

p±i = p±i, m = m + r\ z , (83) 

where the symbols with no hats refer to the lab-frame values. The variables p± and r\ (in addition to the 
asymuthal angle 4>) fully define the lepton momentum, and have an important advantage of very simple 
transformations under boosts along the beamline: p± and ef> are invariant, and rji transforms additively. This 
feature makes this set of kinematic variables extremely convenient, and they are in fact widely used not 
just for the process at hand but throughout hadron collider physics. Their relation to the directly measured 
lab-frame energy E and scattering angle 9 is 

p± = E sin 9 , sinh % = — ^— . (84) 

tan 6/ 
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Fig. 17. Data points: Z boson rapidity distribution measured in pp — > Z — > e + e at the Tevatron by D0. Solid line: NNLO 
theoretical prediction by Anastasiou et.al. 11 From Ref. 13 

Note that the quantity rj z is simply the Z-boson rapidity in the lab frame. In each event, r\ z can be 
determined experimentally: since 771 = —772, we have r\ z = 771 + r\i- On the other hand, 



Thus, measuring the Z rapidity distribution in effect provides the distribution of x values of partons con- 
tributing to Z production, giving an indirect but powerful constraint on pdf's. The distribution measured 
by the D0 collaboration at the Tevatron, along with the NNLO theoretical prediction of Ref., 11 is shown 
in Fig. 17. The agreement between theory and experiment is spectacular. Note that, while the shape of 
the curve depends on the pdf's, the location of the cndpoints is easy to understand: x\ G [xmin, 1] implies 
1*7*1 < 5 logs/Mf = 3.0 at the Tevatron. (At the LHC, ri max = 5.0.) 

The inclusive processes whose cross sections we computed here are sometimes called "W/Z+ > jets" . 
Cross sections with 1 or more jets in the final state are also of interest. In general, cross sections for W/Z+ > 1 
jet, W/Z+ > 2 jets, etc. are easier to compute and compare to experiment than those with exactly specified 
number of jets, due to the ambiguities in counting jets arising from collinear gluon emission (see Section 2.4.3). 
These cross sections depend sensitively on the minimum p± required for the jet(s), since low-pj_ jet rate is 
large due to the collinear singularity in ISR gluon emission. Very roughly speaking, each extra jet in the final 
state with p± <; Mz reduces the cross section by a s ~ 0.1, but the suppression is smaller if p± is lowered. 

3.2.3. W Mass Measurements 

It is easy to measure the Z mass at a hadron collider: One simply needs to select the Drell-Yan events 
(those with with two opposite-charge leptons in the final state), and compute the dilepton invariant mass, 
s i2 = (pi +P2) 2 , in each event. Most Drell-Yan events with large invariant mass come from Z decays, with a 
small contribution from off-shell photon exchanges. While the Z velocity along the beam line varies between 
events, it does not affect S12, which is Lorentz invariant. So, plotting the number of events vs. S12 will produce 
a Breit-Wigner curve, with a peak at Mz and width T z- This method provides a simple and unambiguous 
determination of the mass of any unstable particle, provided that it has at least one decay channel where all 
decay products can be detected, and backgrounds in that channel are manageable. 




(85) 
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Unfortunately, this is not the case for the W. Hadronic W decays produce dijet final states, wh ich ar e 
buried under a large background from pure-QCD events (see Section 3.3). Leptonic decays are clean, but 
involve a weakly-interacting particle, the neutrino, whose momentum cannot be measured. Thus, the simple 
invariant-mass measurement is impossible, and other techniques have to be used. I will outline two useful 
approaches below. 
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Fig. 18. Distributions of the electron pt (left panel) and transverse mass (right panel) in the reaction pp — > W — > eu, measured 
by the D0 collaboration at the Tevatron. From Ref. 14 



The first approach is to look at the p± distribution of the charged lepton. Since (at least to leading order 
in a s ) the produced W is at rest in the parton frame, we have 

M w . - 



Pi 



p± = Esin( 



(86) 



where hats indicate quantities evaluated in the parton frame, and the lepton mass has been set to zero. It 
follows that 

M w 



0<P±< 



(87) 



Moreover, 



da d cos 8 da 



P± 



da 



dp± dcos9 



^ M w ^ 2 _ ^ G?COS( 



It is easy to show that - da ~ does not vanish at 6 = n/2, corresponding to p± — Myy/2. Eqs. (87), (88) then 
imply that the lepton-p^ distribution blows up at its upper boundary, Myy/2, and then abruptly drops to 
zero. This discontinuity is tempered by the W width effects, but a sharp peak at Mw/2, called Jacobean 
peak to reflect its origin in the variable change 9 — > p±, remains. The peak is clearly seen in the distribution 
of electron p± from W decays measured by the D0 collaboration at the Tevatron, shown on the right panel 
of Fig. 18. Measuring the position of this peak determines the W mass. 

A difficulty in using the lepton p± is that the transverse motion of the W , for example arising from recoils 
against an extra jet which can be emitted at NLO in a s , must be carefully taken into account before an 
accurate mass determination can be made. An alternative variable, which does not suffer from this problem, 
is the transverse mass. To define it, note that while neutrino cannot be detected, its momentum in the 
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direction transverse to the beam can be reconstructed, since the sum of all transverse momenta in th e even t 
must be zero by momentum conservation: 

pi = -pi -E pi- (89) 

3 

where the sum is over all jets in the event. (The neutrino momentum along the beamline cannot be recon- 
structed, since the momentum carried by remnants of the colliding protons along this direction cannot be 
measured.) The transverse mass is then defined as 

m\ = (|pl| + |pi|) 2 - (Pl+Pi) 2 . (90) 



Homework Problem 3: Show that < m\ < M w . Hint: compare m\ with the (unobservable) ev invariant 
mass, s ev = (|p"| + |p £ |) 2 - (p" + p £ ) 2 - M w . 

Just like for the lepton p± , the change of variables 6 — > tut introduces a Jacobean factor that blows up 
at the upper boundary vht — M\y, leading to a Jacobean peak around this value. This peak is clearly visible 
in the experimental data, shown in the left panel of Fig. 18. 

A combination of these techniques yields an amazingly accurate measurement of raw at the Tevatron: 
at present, the error is only about 30 MeV, or Am\y/mw ~ 4 x 10~ 4 ! This is comparable to the precision 
achieved at LEP. The lesson is clear: with clever variable choices and thorough understanding of systematic 
issues, hadron colliders may be capable of matching the precision measurements possible in e+e~, even for 
observables which at first glance suffer from the lack of kinematic information. 



3.3. Dijet Production 

The most common processes at hadron colliders are QCD-mcdiated reactions with only strongly interacting 
particles in the final state. Conservation of energy-momentum and fermion number imply that there are no 
2—^1 processes of this kind; the dominant reactions are 2 — > 2 processes, which upon hadronization result 
in two-jet, or dijet, final states. We will briefly consider dijet production in this subsection. 

In general, the only variable needed to describe the kinematics of a 2 — > 2 scattering process in the 
parton frame is the scattering angle 9, or, equivalently, the Mandelstam variable t (for massless particles, 
t = — s(l — cos#)/2). In the case of hadron collisions, there are two additional variables, x\ and x 2 , which 
determine the parton c.o.m. energy as well as the motion of the parton frame with respect to the lab frame. To 
compare with data, it is convenient to change variables from (t, xi,x 2 ) to the directly observable (p± , r/i , rj 2 ) • 
In terms of these variables, the triple-differential cross section has the form 

= 2p ±Xl x 2 J2f il (x 1 )f i2 (x 2 )^(i 1 +i 2 ^ 1 + 2), (91) 



dp±drj 1 dr) 2 dt 

H,l 2 



where 



x\ = cosh ri e n+ , 

X2 = cosh?7_ e~ I,+ , 
t = -2p\ cosh77_e-' 7 - , (92) 



with r]± = (771 ± 772V2. 



Homework Problem 4: Derive the relations (92) and the triple-differential distribution (91). 
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Fig. 19. The Feynman diagram for the process ud — ¥ ud at the tree level. 
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Fig. 20. Jet cross sections in five rapidity regions, as a function of jet pj_, measured by the CDF collaboration at the Tevatron 
and compared with the NLO-QCD calculation. The cross sections for the five regions are scaled by a factor of 10 3 from each 
other for presentation purposes. From Ref. 15 



To compute the dijet rates, we need to know parton-level differential cross sections da/di for each possible 
reaction with two strongly interacting particles in the final state. There is a large number of such reactions. 
Here is a partial list: qq — > qq, where q can be any of the quark flavors; qq — > gg; qq — > q'q' with q' ^ q; 
qq — > qq; qq' — > qq'; gq —¥ gq; qq — > gg; etc. As an illustrative example, consider the subprocess ud — > ud. 
The single Feynman diagram that contributes to this process at leading order in a s is shown in Fig. 19. The 
cross section is 



da 
di 



4na s 



Using s = X1X2S and s + t + u = 0, together with Eqs. (92), it is easy to show that 



da 



1 



-j a — , as p± ->• 0. 
dt p ± 

At the hadron level, the triple-differential cross section (for fixed jet rapidities) behaves as 

d 3<7 3 , / \ , , % 1 

OC p ± ■ fu(CuP±)Jd{CdP±) • — + ... , 

Pi 



(93) 



(94) 



(95) 



dp±driidr]2 

where c u and Cd are constants independent of pj_, and dots indicate contributions from other subprocesses. 
Since the pdf's grow faster than l/iasa:->0, the right-hand side of Eq. (95) grows faster than l/p\, and 
the hadron- level cross section diverges strongly as pi — > 0. This singularity is of a different nature than 
the soft and collinear singularities encountered so far: It is due to an exchange of a massless particle in 
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the t-channel, and is exactly analogous to the Rutherford singularity in the elastic scattering by C oulomb 
interaction in QED. The singularity (and similar divergences that occur in other subprocesses involving t 
or u-channel exchanges) is effectively regulated by the IR divergence of the strong coupling constant: When 
the momentum exchanged in the t channel is of order Aqcd, perturbation theory breaks down and the 
leading-order result is no longer valid. As a result, the total rate for pure-QCD events is simply given by the 
geometric cross section 

47T 

QCD 

This cross section is huge: at current luminosity, there are over 10 7 pure-QCD events per second at the 
Tevatron. This rate is far too high for the experiments to be able to record all events on tape, let alone 
analyze them. However, the vast majority of pure-QCD events involve only hadrons with p± <~ Aqcd, 
whereas events of interest for probing the SM at high energy scales and searching for new physics involve 
much higher transverse momenta. This leads to a strategy called triggering: only events with at least one jet 
above some minimal p± value are recorded and analyzed. (Events with non-strongly-interacting objects, such 
as leptons or photons, are far less numerous than pure-QCD, and most of them are also recorded.) Thus, the 
experimentally relevant cross section is that for producing jets with p± above the threshold required by the 
trigger (or a larger value set by analysis cuts). For example, for p± > 100 GeV, a typical value for the LHC 
analyses, the total jet cross section at the LHC is about 1 mb, about a factor of 20 larger than the inclusive 
Z production cross section. 

Comparing dijet differential cross sections with data provides a valuable test of QCD at high energy scales, 
as well as the most accurate determination of gluon pdf at high x. Within limitations due to theoretical 
and experimental uncertainties in jet definitions, uncalculated higher-order perturbative corrections, pdf 
uncertainties, and other systematic errors, all available data is in agreement with the SM. An example is 
provided by a recent CDF measurement of dijet rates at the Tevatron, see Fig. 20. 

3.4. Hadron Collider Cross Section Summary 

A useful summary of cross sections for the most important SM processes in hadron collisions is shown in 
Fig. 21. The total scattering cross section er t ot ~ 0.1 bn is completely dominated by \ow-p± pure-QCD scat- 
tering, and depends only weakly (logarithmically) on yfs. The inclusive cross section for b quark production 
is 2-3 orders of magnitude lower, reflecting the effective p± cutoff provided by the b mass. The strongest 
electroweak processes, inclusive W and Z production, have cross sections of order 10-100 nb: in other words, 
there is 1 electroweak boson produced per 1-10 million pure-QCD events. Of course, as discussed above, 
the vast majority of pure-QCD events are not recorded on tape, and experimentally relevant QCD events 
are those with at least one jet with p± above a trigger threshold. While specific thresholds depend on the 
experiment and the analysis, the figure shows the total jet cross sections for a few representative values: 
y/s/20, y/s/4, and 100 GeV. As a rough guidance, it is useful to remember that at the Tevatron, the jet rate 
with p± > 100 GeV is roughly comparable to the electroweak boson production rates, while at the LHC it 
is larger by 1-2 orders of magnitude. Going further down the cross section scale, the top pair-production 
cross section is in the few-pb range at the Tevatron, and close to a nb at the LHC. The sharp rise with 
y/s is due to the fact that top pairs are typically produced by partons with x ~ 0.2 at the Tevatron and 
x ~ 0.02 at the LHC; the higher parton density at low x, especially for gluons, increases the cross section. 
At design luminosity, the LHC will produce about 10 7 top pairs a year, allowing for detailed studies of top 
properties (but also providing an important background for many new physics searches). New physics cross 
sections are generally further suppressed: as an example, cross sections for the SM Higgs production, for 
two representative values of the Higgs mass, are shown in the figure. The main challenge in new physics 
searches at hadron colliders is to distinguish the rare new physics events from the much more common SM 
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Fig. 21. Cross sections for a variety of SM reactions, as a function of center-of-mass energies, in pp collisions (for \fs < 4 TeV) 
and pp collisions (for y/~s > 4 TeV). 



processes. To do this, one must identify specific observable features, or signatures, of the new physics events 
that distinguish them from the SM. 



4. Searches for New Physics 

In this Lecture, I will first outline the general strategy for new physics searches in collider experiments. I 
will then consider two specific examples: searches for the Higgs boson, and supersymmetry. 



4.1. General Strategy 

Suppose that you are given a model that extends the SM and predicts new particles at energy scales poten- 
tially within reach of current or near-future colliders. Your task is to plan a search for this model: You need 
to identify its experimental signatures, devise a data analysis strategy in order to maximize sensitivity to new 
physics, and evaluate how much data would be required to convincingly establish, or rule out, the theory. 
(Parenthetically, note that these tasks fall right on the boundary of what theorists and experimentalists do. 
Typically, a first-round study would be done by theorists. If the results look promising, a more detailed 
analysis, taking account of experimental conditions, would be conducted by experimentalists, who will then 
implement the designed search with real data.) Of course, the detailed strategy would depend on the precise 
nature of the predicted new particles and their interactions, as well as on the experiment. However, most 
studies of this kind roughly follow the "twelve-step" blueprint outlined below. 
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(1) Suppose that our model predicts a particle X, and we would like to design a search for this p ar tide . 
Start with the Lagrangian of the model. Identify all couplings of X to the SM. There is no need to derive 
precise Feynman rules yet; simply list the vertices involving X and SM states. 

(2) Given the vertices, identify X production mechanisms at the collider of interest. Draw Feynman diagrams 
with X in the final state. Remember that cross sections generally decrease with the number of interaction 
vertices and with the number of particles in the final state, so focus on drawing diagrams as simple as 
possible. If considering a hadron collider, remember that pdf's drop fast with x, and therefore cross 
sections decrease with the invariant mass of the final state. This implies that the final states with 
highest cross sections will involve as few heavy particles as possible. If the new particle can be singly 
produced, this will be the dominant process; if single production is forbidden by symmetries (e.g. R-parity 
in supersymmetric models), pair-production will dominate. 

(3) Based on the above considerations, identify the most promising production processes for X. You can also 
order diagrams according to the strength of couplings that enter: for example, a s > a w , so, everything 
else being equal, QCD-mediated processes will have higher cross sections than electroweak reactions. 

(4) Carefully derive the Feynman rules for all vertices that enter the promising production processes. 

(5) Compute the total cross section for each of the promising processes. Do not impose cuts, unless necessary 
to avoid singularities; if necessary, impose very loose cuts roughly consistent with experimental conditions 
(if unsure, consult your experimental colleagues). Given the integrated luminosity (actual or expected) 
of your collider, compute the number of events N ev expected in your data set for each of the processes. 
(See Table 1 for luminosity values.) 

(6) If N cv < 10 for a given process, the process is probably uninteresting - do not consider it further. Note 
that iV ev typically depends on a priori unknown model parameters, such as the X mass. As long as there 
is a range of parameter values, consistent with theoretical considerations and with constraints from prior 
experiments, where the iV ov > 10 condition is satisfied, it is worthwhile to continue with the analysis. 

(7) Identify all vertices through which X can decay (both to SM particles and to other exotics). For each 
possible decay channel, check whether it is kinematically allowed. If there are no viable decays, X is a 
stable particle. If X is stable, is it charged under strong and/or electromagnetic interactions? 

(a) If yes, you need to think carefully about the signature of X in the detector. If only EM charge 
is present, X would behave like a heavy muon, and will leave a trace in the muon detector. If X 
is strongly coupled, it will quickly hadronize, and the resulting hadrons with interact in the hadron 
calorimeter, as well as in the tracker if electrically charged. (As an aside, note that there are significant 
cosmological constraints on charged and/or strongly-interacting stable or long-lived massive particles.) 
See K. Zurek's lectures at this school for more details. 

(b) If no, X can at best interact with matter in the detector via weak interactions, or possibly even weaker. 
In this case, X escapes the detector without interacting, like the SM neutrino, and is "invisible". 
Invisible particles can be detected as an apparent missing momentum (transverse momentum in the 
case of hadron collisions) recoiling against visible particles. So, in this case, make sure that there is 
at least one observable "tag" particle in the final state of your process. A jet or a photon, which can 
always be emitted as ISR, can serve as tags. 

(8) If X is unstable (as is usually the case for heavy new particles), list all X decay channels, and compute 
their branching ratios. (Again, focus on the simplest possible decays since those with more interaction 
vertices are suppressed.) If some of the X decay products are themselves unstable, repeat this step 
for those particles, and iterate until you only have stable particles left. List all final states in terms of 
"detector objects" - jets, SM leptons, photons, invisible particles (neutrinos and possibly exotic stable 
invisibles), and possibly exotic charged/colored states. For each final state, compute the probability to 
obtain it from X decay, multiplying the relevant branching ratios. Multiply by N cv to obtain the expected 
signal rate in each channel, iV s i g . 
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(9) For each detector-level final state with a sizable expected signal rate, list all SM processes that c an giv e 
rise to this final state - the Standard Model background. (Again, focus on the simplest possible processes. 
You can use Fig. 21 as a guide to the most common processes and their cross sections, but remember 
that it is incomplete.) 

(10) Compute the cross section for each SM background process, and multiply by the luminosity to get N$m- 

(11) A meaningful search (i.e. one that can lead to a discovery or a publishable constraint on new physics) 
must satisfy three conditions: 

(a) Events: iV s i g > a few. Detector efficiency should be taken into account when evaluating iV s i g ; for 
a first-round, "theory-level" analysis, a rough estimate of the efficiency would suffice. Note that 
A s i g oc Li n t , so this condition can at least in principle always be satisfied by running the experiment 
long enough. 

(b) Statistics: N s i g / \J Nsm > a few. "A few" is typically 3 or 5, depending on how confident one wants 
to be. In case of an actual discovery, 5 would probably be required before most people are con- 
vinced, given the long history of "3-sigma" effects disappearing with increased statistics. In this 
case, TVgig/VAsM oc L\^, so this condition in principle can also be always fulfilled by increasing the 
data sample size, although the improvement is slower than for (a). 

(c) Systematics: N s i g /NsM > (a few) x er sys , where a sys is the fractional systematic uncertainty of the 
SM prediction. In the case of hadron collisions, this uncertainty includes residual scale dependence 
from uncalculated higher-order corrections, pdf uncertainties, as well as experimental errors, e.g. 
uncertainty in the detector efficiency. (Talk to your experimental colleagues to get a reasonable 
estimate for the latter.) Note that N B i g /NsM oc L? nt , so no amount of luminosity will help if this 
condition is not satisfied. Thus, it is crucially important to understand all sources of systematic 
uncertainty, and to eliminate them as far as possible. 

(12) If the three conditions are not satisfied, one can usually improve the situation by imposing selection 
cuts: Look for differences in signal and background kinematic distributions, and devise criteria to se- 
lect only events in the signal-rich and/or background-poor part of the phase space. This improves sig- 
nal/background ratio, making conditions (b) and (c) easier to satisfy. (However, one must be careful since 
(a) can fail if the selection cuts are too tight!) 

Once the three conditions are satisfied, the proposed search strategy can be studied in more detail by 
experimentalists, and, if the conclusions hold up, implemented with real data. You have succeeded! 



4.2. Higgs Searches 

The remainder of this Lecture is a brief overview of some of the most important searches for new physics 
which are the focus of the current experimental program at the Tevatron, as well as the upcoming one at 
the LHC. 

We begin with the search for the Higgs. I will focus on the simplest scenario incorporated in the SM, a 
single Higgs doublet field yielding one physical Higgs boson h. The Higgs potential has two free parameters, 
the mass parameter /i 2 and the quartic coupling A; these can be traded for the Higgs vacuum expectation 
value (vev) v — 246 GeV and the physical Higgs mass m^. The Higgs is coupled to SM fermions via Yukawa 
couplings Ufh'tptp, where yj — \/2mf/v, and to electroweak gauge bosons via the usual gauge vertices. Note 
that tree- level triple-boson couplings of the form hVV exist for massive vector bosons, V = W/Z, but 
not for the photon. The vertices /177 and hgg arise from the one-loop diagrams in Fig. 22; while loop- 
suppressed, these couplings play an important role in Higgs phenomenology. Note that all Higgs couplings 
can be expressed in terms of known SM parameters; the only unknown parameter in this model is rrih- 

Having listed the vertices, we can now identify Higgs production processes at colliders. First, consider 
e+e~ collisions. The simplest reaction would be s-channel resonant production, e+e~ — > h, shown in Fig. 23 
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(a) (b) (c) 



Fig. 22. One-loop Feynman diagrams which induce the Higgs coupling to gluons (a) and photons (b, c). In diagrams (a) 
and (b) contributions of all quarks and leptons in the loop must be summed over, but in practice the top quark contribution 
dominates due to its large Yukawa coupling. 




Fig. 23. Leading-order diagrams for Higgs production in e+e collisions. 

(a) ; however, the cross section of this process is tiny due to suppression by y? e ~ 10~ 12 . The strongest 
production channel is the so-called Bjorken process shown in Fig. 23 (b) , which utilizes the gauge coupling of 
the Higgs. This reaction is kinematically allowed if y/s > rrih+mz- At LEP-2, this implies that the maximum 
Higgs mass that can in principle be probed is about 118 GeV. In practice, the Higgs mass must be a bit lower 
to get a sizable production cross section: LEP-2 expected to produce 10 or more Higgs events for < 115 
GeV. To identify the Higgs signature, we next need to consider Higgs decays. The SM Higgs branching 
ratios, as a function of its mass, are presented in Fig. 24. Throughout the parameter space accessible at 
LEP, the dominant Higgs decay mode is bb, with a sizable (order few-%) contribution from r + r _ . Thus, 
the interesting final states include bb or t + t" pairs in association with lepton or jet pairs from Z decays 1 , 
or missing momentum from invisible Z decays Z — > vv. SM backgrounds for all these channels primarily 
come from WW and ZZ final states, and can be accurately computed. Combining the channels, LEP-2 
experiments place the lower bound on the Higgs mass: 5 

m h > 114.4 GeV, 95% c.l. (97) 

Interestingly, the ALEPH collaboration reported 5 candidate events that may be viewed as a hint for the 
Higgs boson at about 115 GeV; 17 however, low statistics, and lack of confirmation by other LEP experiments, 
malKllmHiigSa^&nch^c^fchs&'^.n colliders are somewhat more complicated. Resonant production by quark- 
antiquark pairs, Fig. 25 (a), is again strongly suppressed by weak Yukawas. The dominant process, at both 
the Tevatron and the LHC, is the one-loop gluon fusion process shown in Fig. 25 (b). However, several 
other processes also lead to useful signatures: the most important ones are the associated production with 
an electroweak gauge boson, Fig. 25 (c), and the vector boson fusion (VBF) reaction, Fig. 25 (d). The cross 
sections for these processes at the Tevatron and the LHC are shown in Fig. 26. The Tevatron should expect 
to collect a sizable sample of Higgs bosons for < 200 GeV, whereas the LHC at design luminosity should 
produce plenty of Higgses throughout the theoretically sensible range, < 1 TeV. The signature of the 



'Such pairs can be identified in the data by the fact that their invariant mass is close to mz- 
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Fig. 24. Branching ratios of the SM Higgs. From the web site. 



Higgs depends crucially on its mass. Roughly speaking, one can divide the parameter space into two regions: 
the low-mass region, below about 140 GeV, and the high-mass region, from 140 GeV to 1 TeV. In the low- 
mass region, the dominant decay mode is bb. The problem is that this mode, combined with gluon fusion 
production, gives a signature of 2 jets, at a rate tiny compared to the SM pure-QCD dijet rate. It is impossible 
to identify the Higgs events on top of this huge background, and one is forced to focus on either subleading 
production processes, or subdominant decay modes. At the Tevatron, the most promising strategy is to focus 
on h — > bb but look for associated production events, Zh and Wh, with leptonic vector boson decays. The SM 
backgrounds in these channels, predominantly Z/W + 2 jets, are much more manageable. At the LHC, where 
the expected Higgs production rate is much higher, one can concentrate instead on a rare but very clean 
decay h — > 77. In the high-mass region, on the other hand, the dominant decay mode is h — > WW (with 
one of the Ws off-shell for rrih < 160 GeV), and h — > ZZ also makes a significant contribution for nih <^ 180 
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Fig. 25. Higgs production processes in hadron collisions: (a) s-channel resonant production (suppressed by small Yukawas); 
(b) gluon fusion; (c) associated production with a vector boson; and (d) vector boson fusion. 
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Fig. 26. SM Higgs production cross section at the Tevatron, in pb (top) and the LHC, in fb (bottom), as a function of the 
Higgs mass. The VBF cross section at the Tevatron is very small, and is not shown here. The Tevatron plot is from the web 
site, 19 and the LHC plot is from Ref. 20 



GeV. If the vector boson decays leptonically, these modes do not suffer from large QCD backgrounds, and 
provide a clean signature. The "golden" h — > ZZ — > A£ mode is especially advantageous since it contains no 
invisible particles, so the Higgs boson shows up as a clear peak in the 4-lepton invariant mass distribution. 
Combined with the strong gluon fusion production channel, the vector boson decays make the search in the 
high-mass region relatively straightforward, giving the Tevatron and especially the LHC sensitivity to high 
Higgs masses in spite of the decreasing production cross section. 
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Fig. 27. Tevatron bounds on the SM Higgs production cross section, as a function of the Higgs mass. The bound combines 
searches for a variety of signatues (see text for details), at both CDF and D0. From Ref. 18 

The latest published results of the SM Higgs searches at the Tevatron are summarized in Fig. 27. To 
present the negative results of the searches, the experiments assume that the Higgs production cross section 
in each channel is a free parameter, and put an upper bound on this cross section. (If several production 
channels are used in a search, the ratios of their cross sections are assumed to be the same as in the SM.) This 
bound, expressed in units of the SM production cross section for each Higgs mass, is plotted on the figure. 
The SM Higgs can be said to be ruled out if the cross section bound in these units is below 1. Currently, the 
Tevatron experiments only exclude the SM Higgs at the 95% c.l. in a narrow mass window just above the 
h — > WW decay threshold, between 163 and 166 GeV. As the integrated luminosity increases, the Tevatron 
is expected to expand this window, unless a deviation from the SM is seen. 

LHC experiments are expected to be sensitive to the entire theoretically consistent range of the Higgs 
masses at the 5-sigma level with about 30 fb _1 of data, see Fig. 28. This corresponds to 3 years of running at 
design luminosity, although of course depending on the Higgs mass such sensitivity may be achieved much 
sooner. Note that, just as at the Tevatron, the LHC sensitivity is not a monotonic function of the mass: in 
fact, a 300-GeV Higgs is substantially easier to discover, via the h — » ZZ — s- 4£ mode, than a 115-GeV Higgs 
where one must rely on the rare decay h 77. 

4.3. Searches for Super symmetry 

As every TASI student knows, the discovery of the Higgs would mark the beginning, not the end, of the 
hunt for new physics at the TeV scale. Strong theoretical arguments indicate that the SM Higgs cannot 
live in isolation, requiring other new particles to stabilize its mass and complete the electroweak symmetry 
breaking sector. The nature of these new particles can at present only be guessed at. One popular guess is 
that physics at the TeV scale is super symmetric: for each particle in the SM, a new "superpartner" particle 
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Fig. 28. Expected sensitivity of the CMS detector at the LHC to the SM Higgs, as a function of the Higgs mass. From Ref. 1 



appears at this scale, with the same quantum numbers but spin different by h/2. Appearing in loops, these 
particles would cancel the notorious quadratic divergence of the Higgs mass parameter present in the SM, 
restoring naturalness to the Higgs sector. I will not have the opportunity for a thorough review of the model 
here; a reader in need of such a review is referred to the set of lectures by P. Fox at this school, as well to the 
excellent review article by Martin 21 and the comprehensive book by Baer and Tata. 22 Instead, I will simply 
very briefly mention the features of the model crucial for collider searches. 



Table 3. The undiscovered particles in the Minimal Supersymmetric Stan- 
dard Model. From Ref. 21 
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The particle content at the TeV scale predicted by the Minimal Supersymmetric Standard Model ( MSSM ) 
is summarized in Table 3. There arc 34 new particles waiting to be discovered! Notice the appearance of 3 
additional Higgs boson states not present in the SM; this is due to the presence of a second Higgs doublet 
field, required in the MSSM for theoretical consistency. All other listed particles are simply superpartners 
of known SM states. If supersymmetry were an exact symmetry of nature, all superpartners would have 
exactly the same mass as their SM counterparts; of course, this possibility is already completely ruled 
out by data. Thus, SUSY must be broken, lifting the superpartner masses into the TeV domain. (Larger 
masses would render SUSY irrelevant for solving the hierarchy problem.) Many models have been proposed 
to describe SUSY breaking (some promising directions were described by P. Meade and D. Shih at this 
school). However, despite over two decades of intense theoretical work, no completely compelling model has 
yet emerged. In the absence of clear theoretical predictions, a conservative approach is to simply consider a 
completely general SUSY breaking Lagrangian, imposing only the restriction that SUSY breaking must be 
"soft" , to not interfere with the SUSY solution to the gauge hierarchy problem. Such general soft Lagrangian 
contains about 100 free parameters, including all superpartner masses, as well as certain trilinear couplings 
among scalar superpartners (so-called "A-terms"). The superabundance of free parameters is a curse on the 
study of SUSY signatures at colliders: Search strategies should ideally be designed that cover as much of 
the huge parameter space as possible, and, equally daunting, the results of each search should be presented 
as exclusion regions in this space. In practice, this challenge has not been met: most searches are planned 
and executed in terms of simplified models that impose relations among parameters, with the 5-parameter 
"minimal supergravity" , or mSUGRA, model being the most common choice. 

In spite of the many free parameters, it is possible to identify generic qualitative features of collider events 
with superparticles. To avoid rapid proton decay, the MSSM postulates an additional discrete symmetry, 
the R-parity. This is a Zi symmetry, i.e. any physical particle is either odd or even under it. Every term in 
the Lagrangian must be even, so that an even number of R-odd particles must appear in each interaction 
vertex predicted by the theory. All known SM particles, as well as all Higgs scalars, are postulated to be 
R-even. The R-parity is embedded in the SUSY algebra in such a way that superpartners must have opposite 
R-charges; in other words, all superpartners of the SM particles are odd. This has two consequences, both 
crucial for collider searches: 

• Superpartners must be pair-produced in collisions of SM particles. 

• The lightest superpartner (LSP) must be stable, and any other superpartner must decay to the LSP plus 
(any number of) R-even particles, either SM or extra Higgses. 

Strong cosmological bounds on electrically charged and colored stable particles indicate that, if the MSSM 
is indeed realized in nature, the LSP should be one of the uncharged, uncolored particles of the model. There 
are three possibilities: the lightest neutralino N\, the lightest sneutrino z>, or the gravitino G. Of the three, the 
neutralino is an especially interesting possibility, since it can easily have the correct relic abundance to explain 
the observed dark matter, without conflict with direct dark matter search experiments. (The neutralino also 
naturally emerges as the LSP in several popular models of SUSY breaking, including mSUGRA.) All these 
particles interact with ordinary matter too weakly to leave a trace in collider detectors, and so they would 
appear invisible, like the SM neutrino. Thus, the generic signature of SUSY at colliders is the excess of 
events with "missing momentum" , accompanied by the SM particles from decays of the originally produced 
superpartners to the LSP j . This qualitative prediction depends only on two basic assumptions -exactly 
conserved R-parity and neutral LSP - and is independent of the precise values of the model parameters. 
Thus, quite generally, an inclusive search for events with missing energy in association with jets and/or 



•>Onc should keep in mind that if the LSP is the gravitino, the next-to-lightcst superpartner (nLSP) may be sufficiently long- 
lived to decay outside the detector. If the nLSP is charged or colored, in this scenario one would look for exotic heavy particles 
leaving traces in the detector, instead of the usual missing energy signature. 
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charged leptons is the best model-independent search strategy. 




Fig. 29. Leading-order diagrams for the MSSM scalar electron production in e+e collisions. 




Fig. 30. LEP-2 bounds on the selectron mass (left) and the sbottom mass (right), plotted vs. the neutralino LSP mass. From 
the LEP2 Joint SUSY Working Group web site. 23 

In e + e~ collisions, sleptons, squarks and charginos can be produced via s-channel 7* jZ exchanges. In 
addition, selectrons, charginos and netralinos can be produced via i-channel exchange of neutralinos, sneu- 
trinos and selectrons, respectively. As an example, consider the selectron search. The Feynman diagrams 
for selectron production are shown in Fig. 29. Produced selectrons can decay directly to electrons and the 
lightest neutralino, assumed to be the LSP. This leads to the e + e~+missing energy final state, which can 
be easily identified at LEP. The SM backgrounds in this channel are not large, and the search is sensitive 
as long as the selectron production cross section is significant. Kinematically, selectron pair production is 
possible if 2m(e) < y/s; this means that LEP-2 could in principle search for selectrons up to about 104 GeV. 
In practice, as usual, the bound is somewhat lower: 

m(e)>100GeV. (98) 

However, there is a small caveat. The search requires that the electrons and positrons from selectron decays 
be identified in the detector. If the selectron mass is close to the LSP mass, these particles are very soft 
and cannot be registered, so the search sensitivity is lost for m(e) — m(JVi) < 1 GeV. The region of the 
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m{e)/m{Nx) plane excluded by the LEP-2 experiments is shown in Fig. 30. Note that this plot is quite robust 
with respect to variations of the other MSSM parameters, since the bound is primarily determined by the 
kinematics of selectron production and decay. The searches for all other electrically charged superpartners 
(sleptons, squarks, and charginos) proceed in a similar way, and all yield bounds close to 100 GeV. For 
example, the LEP-2 bound on the sbottom mass (along with, for comparison, the CDF bound on the same 
parameter) is shown in Fig. 30. 




Fig. 31. Yellow (shaded) bands: Theoretical cross section for total squark/gluino production at the Tevatron, as a function of 
the gluino mass, assuming squark mass of 460 GeV (left panel), or equal squark and gluino masses (right panel). Triangles: the 
upper bound on the cross section observed by the CDF collaboration at the Tevatron. The bound assumes that squark/gluino 
decays follow the pattern predicted by the mSUGRA model. From the web page 24 (see also Ref. 25 ). 



At a hadron collider, the dominant superpartner production process is the production of squarks and 
gluinos via strong interactions: qq* , qq, q*q* , qg, q*g, and gg final states are all possible. The cross sections for 
these reactions depend only on squark and gluino masses, and thus can be robustly predicted. The combined 
cross section for all squark/gluino production channels at the Tevatron is shown by the yellow bands in 
Fig. 31. A large sample of squark and gluinos could in principle be produced at the Tevatron: for example, for 
m(g) — 200 GeV, the current data sample would contain of order 10 5 such particles. As usual, the produced 
squarks and gluinos would decay promptly; the simplest decays are q — > qN and g — > qqN , resulting in 
2-jet, 3-jet, and 4-jet final states associated with missing transverse momentum. (In this context, missing 
transverse momentum is often referred to, somewhat misleadingly, as "missing transverse energy" , or MET. I 
will reluctantly follow this practice.) The SM background to this signature primarily comes from the reaction 
Z + n jets, with invisible Z decays, Z — ¥ vv. In addition, instrumental backgrounds, from pure-QCD multi-jet 
events where apparent MET results from mis-measurement of one or more of the jets, is also important due 
to very large pure-QCD event rates. All backgrounds, however, fall much faster than the expected SUSY 
signal with increased MET, as shown in Fig. 32. (The amount of MET in SUSY events is controlled by the 
squark-neutralino or gluino-neutralino mass difference, of order few hundred GeV; the backgrounds peak at 
0.) Imposing a lower cut on MET in the analysis improves sensitivity to the SUSY signal. Optimizing this 
cut, the CDF collaboration 25 obtained the bounds on the squark/gluino production cross section shown in 
Fig. 31. This was translated into a bound on squark and gluino masses, Fig. 33. Roughly speaking, gluinos 
below 300 GeV, and squarks below 350 GeV, are excluded in this scenario. However, note that while the 
squark and gluino production cross sections are rather robust, their decay chains, and the kinematics of 
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Fig. 32. The MET distribution of events with > 2 jets and MET. Solid line: SM prediction. Dashed line: SM+SUSY (mSUGRA) 
prediction for 349 GeV gluino and 385 GeV squark. Data points: measurement by the CDF collaboration at the Tevatron. From 
the web page 24 (see also Ref. 25 ). 




the decay products, depend on numerous MSSM parameters beyond their masses: most notably, the LSP 
mass is important, since it controls both the amount of MET and the typical jet momenta in SUSY events. 
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The CDF analysis 25 assumed mSUGRA relations between the MSSM parameters; this implies, amongrrrticr 
things, the ratio of the LSP and gluino masses of about 1:7. The cuts were optimized with this assumption. 
It should be noted that, if the MSSM is correct but the mSUGRA is not, such optimization could result in a 
loss of sensitivity: in other words, it is possible that the data contains a hint for SUSY that would be missed 
by this analysis. 26 Such considerations should be taken into account when designing future SUSY searches. 
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Fig. 34. Discovery reach for SUSY projected by the ATLAS collaboration at the LHC. The mSUGRA model is assumed, and 
the reach is shown in terms of the mo and mi/2 parameters, with the other parameters fixed: tan = 10, Ao = 0, fJ, > 0. The 
reach is shown for n-jet+MET signatures. Integrated luminosity of 200 pb — 1 (1/50 of a year at design luminosity) is assumed. 
From Ref. 27 



The same search can be repeated at the LHC. It will benefit from a much larger production cross 
section: inclusive SUSY cross sections in the tens of pb range (or 100,000's events/year at design luminosity) 
are possible in parameter regions not already ruled out by the Tevatron. The sensitivity of the search 
projected by the ATLAS collaboration is shown in Fig. 34. Even with only 200 pb -1 of data, the LHC will 
dramatically increase the reach of the Tevatron, so very early discoveries are possible. (It should be kept in 
mind, however, that unambiguous measurement of MET requires that all potential instrumental sources of 
MET, e.g. cracks in the detector, mis-calibrated calorimeter cells, etc., must be identified first. This task may 
require substantial additional amount of data.) Eventually, the LHC is expected to cover the entire range 
where SUSY provides an attractive solution to the hierarchy problem: If it is not discovered, then either 
electroweak symmetry breaking is severely fine-tuned, or some other mechanism for stabilizing the Higgs 
mass must be found. 



5. Further Reading 

Collider physics is a vast subject, and in four lectures it is only possible to touch the tip of the iceberg. 
Luckily, a variety of excellent resources is available to students wishing to delve deeper into the subject. 
The popular Quantum Field Theory textbook by Peskin and Shroeder, 1 which I referred to throughout 
these lectures, contains an insightful discussion of the key theoretical concepts crucial to the field. Sadly, 
some of these concepts are often not given sufficient emphasis in a typical one- year QFT course: Infrared 
divergences (chapter 6) and evolution equations (section 17.5) are two examples that come to mind. I tried 
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to give a flavor of these topics in these lectures, but interested students are strongly encouraged t o stud y 
them further. At a slightly more advanced level, the text by Ellis, Stirling and Webber 28 contains a wealth 
of material on applications of QCD to e + e~ and hadron collisions, as well as on electroweak processes and 
Higgs searches at hadron colliders. The text by Barger and Philips 29 covers similar subjects, in addition to 
collider phenomenology of a few extensions of the SM. For readers primarily interested in physics beyond 
the SM, the book by Baer and Tata 22 provides a comprehensive discussion of SUSY searches. Another useful 
source is the recent review article by Feng et. al. 30 Pedagogical introduction to collider phenomenology 
of non-SUSY extensions of the SM is harder to find, although there are useful review articles on at least 
some of the popular models (e.g. models of dynamical electroweak symmetry breaking, 31 Little Higgs, 32 and 
Randall-Sundrum and Higgsless models 33 ). In addition to theoretical work, anyone interested in the subject 
ought to pay attention to results coming out of the currently running Tevatron experiments. Both CDF 
and D0 collaborations maintain web pages 34 ' 35 that provide easy and convenient access to both published 
and preliminary public results, grouped by subject. For the LHC, ATLAS and CMS collaborations have 
published Technical Design Reports (TDRs), 36,37 which describe detailed physics studies performed by the 
collaborations to assess their potential for a large number of measurements and searches. For a theory 
student, studying these Reports is an excellent preparation for understanding the LHC data. 

Note Added: In the Fall semester of 2009, I taught a class on Collider Physics at Cornell, which covered 
many of the topics discussed here in much more detail. All (handwritten) lecture notes can be downloaded 
from the course web page. 38 
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